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Abstract 

We discuss the i^-matrices associated to the i?-matrix of a general A^-state vertex model 
whose statistical configurations encode A^ — 1 U{1) symmetries. The factorization condition is 
shown for arbitrary weights being based only on the unitarity property and the Yang-Baxter 
relation satisfied by the i?-matrix. Focusing on the = 3 case we are able to conjecture the 
structure of some relevant twisted monodromy matrix elements for general weights. We apply 
this result providing the algebraic expressions of the domain wall partition functions built 
up in terms of the creation and annihilation monodromy fields. For = 3 we also exhibit a 
i?-matrix whose weights lie on a del Pezzo surface and have a rather general structure. 
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1 Introduction 



The R-matrix plays a fundamental role in the construction of two-dimensional integrable sys- 
tems of statistical mechanics. This operator represented here by Rabija, ^b) acts on the tensor 
product of two A^- dimensional vectors spaces K ® H depending on the complex parameters 
and ^t). The i?-matrix is required to satisfy the Yang-Baxter equation [1], 

i?12(6,6)i?13(ei,e3)i?23(e2,6) = i?23(6,6)i?13(ei,6)i?12(ei,6)- (1) 

The inverse of the i?-matrix can be assured by imposing the unitarity condition, 

i?i2(ei,e2)i?2i(e2,ei)=Xi®x2. (2) 

where is the N x N identity matrix in V^. 

It turns out that the tensor products of the i?-matrices called monodromy operators are central 
objects in the theory of integrable systems [2l|3]. In recent years, it has been realized that such 
monodromy matrix can be decomposed in a suitable way by means of auxiliary operators that 
have been denominated F- matrices [1]. This concept was originally introduced for the six- vertex 
model motivated by the notion of twist deformations of quantum groups [5j. Lets us denote 
by r\'^\{C,i, . . . the product of i?-matrices associated to an arbitrary permutation a of the 
symmetry group S^. The factorization condition for the invertible F-matrices defined by any 
element a{l, ■ ■ ■ ,L) = {o"(l), . . . , cr(L)} G Sl reads as [HE], 

-^cr(l)...cr(L)(^<7(l), • • • ,^(t{L))-Ri!^^l(6, ■ ■ ■ = Fi i{^i, . . . (3) 

where the F-matrices Fi ^(^i, • • • , ^l) act on the tensor product spaces Vi • • ■ (g) Vl- 

The F-matrices can be used as a natural basis to transform the monodromy matrix in a way 
that it becomes totally symmetric with respect to a general permutation of the indices 1 . . . L. 
This similarity transformation for the six-vertex model permits the development of an alternative 
approach [6j to deal with the combinatorial problem underlying the general theory of the scalar 
product of Bethe states [71IH] and the respective computation of domain wall partition functions 
[9]. In some respect this method paved the way for further progress on the formulation of the 

1 



correlation functions for the spin- 1/2 Heisenberg chain [T0HT2] . It also prompted the search 
for explicit forms of F-matrices associated to other integrable vertex models such as for certain 
generalizations of the six- vertex model p^lT^ as well as for multi-state vertex models whose 
weights are based on the SLlnlm) superalgebra [TSHTT] . We also remark that a diagrammatic 
interpretation of the factorization equations for the symmetric six- vertex model has been discussed 
in [18]. 

Recently, it has been argued that the existence of the F-matrices for an arbitrary six-vertex 
model can be pursued without the need of using any explicit weights parameterization [19]. The 
structure of the F-matrices depends basically on the statistical configurations encoded in the R- 
matrix and the verification of the factorization condition can be done by using the algebraic 
weight constraints derived from the Yang-Baxter and unitarity ([2]) relations. We think that this 
point of view of considering the formulation of F-matrices should not be particular to the six- vertex 
model. In this paper we show that this framework can indeed be generalized to tackle integrable 
N-state vertex models that are invariant by — 1 U{1) symmetries. Recall that for N = 2 one 
obtains the standard asymmetric six-vertex model. We apply the aforementioned construction to 
the next simplest case = 3, presenting the algebraic expressions of relevant monodromy matrix 
elements in the F-basis and the corresponding domain wall partition functions. 

This paper is organized as follows. In the next Section we define the U{1)^'^~^^ invariant 
vertex models and write the algebraic relations f|T][2|) for the Boltzmann weights. These explicit 
relations are required to carry out simplifications independent of parameterizations. We motivate 
our approach by exhibiting a solution of the Yang-Baxter equation for = 3 which contains a 
number of free parameters. In Section [3] we discuss a procedure to build up the F-matrices for an 
arbitrary U{1)^'^^^^ vertex model. It combines, in an effective way, past formulations of the F- 
matrices for specific weights with a recent construction devised for the six- vertex model |19j . 
We use the F-basis in Section H] to provide the expressions of certain relevant monodromy matrix 
elements for the = 3 vertex model with general weights. In Sections [5] and [6] we apply these 
results to exhibit the domain wall partition functions associated to products of creation and 
annihilation fields. Our conclusions are presented in Section [71 In Appendices A-C we summarize 
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technical details helpful for the understanding of the main text. 



2 The U{1)^-^ vertex model 

Consider a vertex model whose statistical configurations on both horizontal and vertical links 
of a square L x L lattice take values on possible states. As usual the corresponding row-to-row 
transfer matrix can be written as the trace over an auxiliary space Aa of the monodromy operator 
Ta,i...Lif^)- This matrix is constructed by the following ordered product of /^-matrices, 

7^,1. = RaLifJ-, ^L)-Ra(L-l)(/^, ^L-l) • • • -Ral(/^, 6)- 

From the local Yang-Baxter equation ([T]) it follows that the monodromy matrix satisfies the 
following global intertwining relations called Yang-Baxter algebra, 

In this paper we shall be considering a particular family of A^-state vertex models whose 
statistical configurations are invariant by A^ — 1 f/(l) symmetries. We shall denote the local 
generators of such f/(l) symmetries by S^^'^\ i = 1, . . . , N — 1. This means that the corresponding 
i?-matrix is constrained by the commutation relations, 

[^i2(ei,6),5f''^®X2 + Xi®5^*)] = , t = l,...,N-l. (5) 

In terms of the Wely N x N matrices, e^"^^ G Vj, the expressions for the A^ — 1 azimuthal spin 
operators Sj^'^^ are, 

gi.,^> ^ _ ^+m+D) ^ , = (6) 

Taking into account the property outlined in Eqs.( l5|6|) one finds that the _R-matrix -Ri2(^i,^2) 
has N{2N — 1) non-vanishing weights. For N = 2, the operator Sj^'^^ reduces to the third 
component of the spin 1/2 Pauli matrices and the corresponding statistical system is the fully 
asymmetrical six-vertex model. The possible statistical configurations for general A^ are given in 
terms of three distinct classes of weights denoted here by 04(^1,^2)5 bij{(,i,(,2) and Cjj(,^i,^2)- In 
Figure [T] we show the respective vertex configurations on the square lattice. 
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Figure 1: Elementary configuration of Boltzmann weiglits. 

From Figure [T] we see tliat tlie expression of the associated i?-matrix in terms of of the Weyl 
matrices is given by, 

N N 



^12(6,6) 



i=l 



(ii) 
2 



(7) 



We recall that the integrable vertex model given by Eq.Q with parameterized weights has 
been considered in the literature for some time [20ll2T] . We stress however that the main results of 
this work will be established without the need of any specific parameterization of the Boltzmann 
weights. We shall rely solely on the algebraic relations for the weights aj(,^i,^2), and 
Cjj('Ci) ^2) coming from the Yang-Baxter and unitary properties. For that reason we need to quote 
them here explicitly. By substituting the expression for the i?- matrix ((71) into Eq.(l2]) we find that 
the Boltzmann weights are required to satisfy the three distinct types of relations, 



1 , I 



l,...Ar 



^^j = l,...,iV 



&ij(6,6)cji(6,6) + Cij(^i,6)&»j(6,6) = , ii- 3 = \ 



N. 



(8) 
(9) 

(10) 



The Yang-Baxter equation ([T]) generates extra constraints on the Boltzmann weights depend- 
ing now on three independent rapidities. By substituting Eq.([7]) into Eq.(IT]) one finds that the 
aforementioned functional relations are much more involved. They can however be written in a 
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compact way by the following expressions, 

&^,(6,6)&ife(6,e3) = &,ik(6,6)&i,(6,e3) i 7^ A: (12) 

C3)&ifc(6, 6) = C3)&ifc(6, 6) iy^J 7^k (13) 

^'ij(6, 6)«i(Ci, 6)Qj(6, 6) + 9i(6, 6)cij(6, C3)biji^2, 6) 

= ai(6, 6)&ii(ei, e3)cij(6, 6) i 7^ J (14) 
bijiCi, 6)ai(Ci, 6)cji(6, 6) + CijiCi, 6)cji(6, C3)bij{C2, 6) 

= ai(6, 6)&y (6, 6)cji(6, 6) « 7^ i (15) 
&ii(Ci, 6)Qj(6, C3)&ij(6, C3) + Cij(^i, 6)«-K0. C3)Qi(6, 6) 

= ai(6, 6)cij(Ci, 6)ai(6, 6) « 7^ j (16) 

^ii(6, 6)Qi(6, C3)&ii(6, C3) + Qi(6, 6)aj(6, 6)cij(6, 6) 

= %(6, 6)cij(6, C3)aj(6, 6) ^ 7^ J (17) 
Qj(Ci, 6)ai(6, C3)&ji(6, 6) + ^'ji(Ci, 6)Qj(6, C3)cji(6, 6) 

= Q,(6,e2)&,.(6, 6)0.(6, 6) ^ 7^ J (18) 

Qj(6, 6)%(6, 6)^u(6, 6) + ^ij(Ci, 6)Qj(Ci, 6)9^(6, 6) 

= Q,(6, 6)&^,•(6, 6)%(6, 6) i 7^ j (19) 

Cij(6, 6)cjfc(6, 6)%(6, 6) + %(6, 6)cifc(6, 6)cii(6, 6) 

= Qfe(ei, 6)i'i,(a, e3)c,ik(6, 6) « 7^ ^ (20) 

Cifcj(6, 6)Qit(6, 6)^'jfc(6, 6) + &ifc(6, 6)cij(6, 6)9^(6, 6) 

= Cij(6, C2)&jfc(6, 6)cifc(6, 6) i^jj^k (21) 

^'y (6, 6)Qfe(Cl, C3)&ji(6, 6) + Qj(Cl, 6)9fc(6, C3)Cy(6, 6) 

= &fcj(6, 6)cifc(Ci, ^s)bjk{^2, 6) + Cjfc(6, 6)cii(6, 6)cjfe(6, 6) ^T^JT^k (22) 
where the indices i,j,k — l,...,N. 
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We would like to conclude this section by remarking that the Yang-Baxter Eqs. (llH l 22|) hide 
integrable models whose weights structure are indeed rather general. Recall here that the case 
N = 2, the asymmetric six- vertex model, has already been detailed in the literature [ni22]. 
Therefore, we shall concentrate on the next simplest system which is the = 3 fifteen-vertex 
model. The key ingredient to solving integrable models whilst keeping their weights as arbitrary 
as possible is to uncover the main algebraic varieties constraining the respective Boltzmann weights 
m . In order to tackle this problem for the case = 3 we have adapted a method, first developed 
in [23], which handles a large number of functional equations associated with three-state vertex 
models. The technical details of this analysis have been summarized in Appendix A and in what 
follows we will present only the main results. It turns out that the underlying algebraic variety of 
one possible solution is a homogeneous hypersurface given by the equation. 



where Ai,A2 are free constants while a(^i), b{^i), c(^j) and c(^j) are arbitrary variables 

depending on the spectral parameters. 

One possible way to parameterize the hypersurface (I23p is first to consider its intersection with 
the hyperplane c{^i) —c{^i) = 0. As a result we obtain an algebraic variety in the class of the cubic 
del Pezzo surfaces which can be parameterized in terms of rational functions [24J. Following an 
algorithm devised in we conclude the rational map is attained by just solving Eq. fl23l) for the 
linear variable in terms of the remaining parameters ct(^j),&(^j) and c(^j). Clearly, the same 
type of procedure also works for the general manifold (l23ll and therefore the i?-matrix depends at 
least on the three free variables and c(^j). Taking into account the results of Appendix 

A we find that the i?-matrix elements are given by. 




(23) 



Cl2(^l,6) 



c(ei) [a(6) - Ai6(6)] [(A1A2 - l)a(6) - A^^)] 



(24) 



^12(^1, 6) 

Cl2(6,6) 




+ AiA2)c(6)c(ei) [a(6)K6) - a(^i)K6)] 

[-a(6) + AM^2)r' [(A1A2 - l)a(6) - Ai&(e2)] 
[-a(ei) + AiK^O] [(A1A2 - l)a(ei) - AMi)] 



1 



X 



(25) 
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fei3(6,6) 

Cl2(6,6) 



[(Ai A2 - l)a(6) - Ai6(6)] [(Ai A2 - l)a(6) - Ai6(6)] 



(26) 



fe2l(G,6) 
Cl2(^l,6) 



a(e2)fe(6) - a(6)K6) 
c(6)c(6) 



a2(6,6) ^ c(6) [(A1A2 - l)a(6)a(6) - Af&(6)(A2a(ei) - K^i))] .28' 
ci2(ei, 6) c(6) [a(a) - Ai6(a)] [-a(6) + ^i^2a{^i) - A^i)] ' 

&23(ei,6) ^ (A1A2 - i)/ii(6)/^2(ei) [«(6)Kei) - a(G)&(6)] . . 

ci2(6,6) Mei)c(6)c(6) [(AiA2-l)a(6)-Ai6(ei)] ^ ' 

ci3(ei,e2) ^ c(6)c(6)/^i(ei) [a(6) - AiK^i)] 

ci2(6,6) c(ei)c(ei)/ii(6) [a(6)-Ai6(6)] ^ ■ 

C2l(6,6) _ C(6)c(ei) C23(6,6) _ c(6)/^l(ei) C3i(ei,6) _ c(6)/^2(ei) .3.- 

012(6,6) c(6)c(6)' ci2(ei,6) c(6)/ii(6)' ci2(6,6) c(6)/^2(6) ^ ■ 

&3i(6,6) ^ ^2c(6)c(6) H^2)m) - «(6)K6)] . . 

ci2(ei,6) /^i(6)c(6)/^2(6) [a(6) - Ai6(e2)] ^ 

032(6.6) ^ /^2(6) [(A1A2 - l)a(6) - Ai6(6)] 

ci2(6, 6) /^2(6) [(A1A2 - i)a(6) - Ai6(6)] ^ ■ 

&32(6,6) _ Af(5ic(6)c(6)c(6) [a(6)K6) - «(6)K6)] 

ci2(6, 6) /^i(6)/^2(6) [(Ai A2 - l)a(6) - Ai6(6)] [a(6) - Ai5(6)] [a(6) - Ai6(6)] 



«3(6.6) _ /^i(6)c(6)c(6)/^2(6) [(A1A2 - i)a(6)a(6) - Af&(6)(A2a(6) - K6))] 
ci2(6,6) /^i(6)c(6)c(6)/i2(6) [(A1A2 - l)a(6) - Ai6(6)] [a(6) - Ai6(6)] 



where 012(6,6) is an overall normalization, and 62 are extra free constants while /ii(6) and 
^2(6) are free variables. The latter freedom is related to the fact that the Yang-Baxter equation 
is preserved under local transformations associated to the two U{1) symmetries. 

The structure of the i?-matrix is certainly more general than that of the the standard U{1) ® 
U{1) invariant vertex models containing only one free spectral parameter [201121] • This rather 
particular parameterization for the completely free variables 0(6)5^6) and c(6) is discussed 
at the end of Appendix A. This fact highlights the importance to search for explicit results in 
integrable models that are independent of specific parameterization of the Boltzmann weights. 
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3 Factorized F-matrices 



We begin by providing some basic definitions and notation that are going to be used in the 
text. We first recall that an element a of the symmetric group Sl can be generated in terms of 
adjacent permutations. In other words one can write, 

<7 = CrQp(ap + l) • • • 0"Q2(a2 + l)'^Ql(Ql + l)5 (36) 

where aa{a+i) denotes the permutation of the indices a and i.e. aa{a+i) (1, . . . , a, a + 1, . . . , L) 
(1, . . . , a + 1, a, . . . , L). We shall refer to Eq. (136|) as the minimum decomposition of a. 

In order to avoid cumbersome notation by the presence of the inhomogeneities in the F- 
matrices we shall omit them when this does not generates confusion. In general, we shall identify a 
given general element Xa,a(i),...,a(L)(/i|^<7{i), • • • ,^<7{l)) e End(A(S)V;(i)®- • -OKcl)) to Xa,^(l...L)(/i)• 
Taking this notation into account the factorization condition (|3]) is rewritten as, 

F^(i-L)R''A = Fi...L. (37) 

The tensor product of i?- matrices R^^i entering Eq. (j3711 is defined in terms of product of 
auxiliary operators by the expression, 

^11 = P^IR^:!^- (38) 

Through considering the minimum decomposition of a in terms of adjacent permutations, the 
auxiliary operators Pi'^\ and B!f ^ ^ can be written in the following way, 

1...L ~ 1...L • • • 1...L 1...L /oq\ 

pi^^"^} _ p{'^"l("l+l)} p{'^"2(<»2 + l)} p{<^«p(Qp + l)} 

^l...L — ^l...L ^l...L ••• .L ' 

where -P^..^^"^^'^ = Pa(a+i) is the standard permutator while Ri,,]^"^^^^ = Pa{a+i)Ra{a+i)- 

We now list a number of properties that are necessary in our analysis of the factorization 
condition (157|) . We start by mentioning that for any given operator Xi ^ £ End(Vi ® ■ ■ ■ Vl) 
we have the following useful relation, 

^.(i...l) = pSXi...lp£'^. (40) 



Next, given the decomposition laws of the auxihary operators (139|) one can show that for 
general {cr, r} G Sl the tensor product of permuted i?-matrices obeys the following decomposition 
identity, 

= (41) 

In order to verify Eq.( HT]) we consider the expression R\^^1 and perform the following opera- 
tions, 

pW _ pW} A{(ar)-i} _ pW pM Mr-'} p{<x-i} 

^^2^^ ~ ^l-.-L ■^l...L — ^l...L fl...L-'^l...L ^ -"-l...L 
apply Eq.|[38t apply Eq.|[381| 

= ^ = (42) 

apply Eq.gnil apply Eg. 13811 

We end by mentioning another standard identity regarding the action of the operator rI'^^ on 
the monodromy matrix, namely 

i?Jir.,i..x(/^) = r.,.(i...L)(/i)i?Jl. (43) 

3.1 The L = 2 example 

We shall start the explicit construction of the F-matrices for arbitrary weights 0-1(^1,^2); 
^2) and Cjj(^i, ^2)- Besides the factorization relation (1S7|) we also require that the F-matrices 
are lower triangular and invertible jlj. In this sense it is instructive to begin with the simplest 
case L = 2. In this situation the only non-trivial permutation is the adjacent permutation au 
and Eq. (l37I) becomes, 

F21RU = Fu. (44) 
We shall now explicitly show that the solution the expression is given by. 



Fi2 = AfuJ'u (45) 

where A/12 is the diagonal matrix, 

N N 

=Xi®X2-Y. ® ^ + E V«^(^i>6)ef^ ® ef ) , (46) 



i=l i=l 
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and Tyi is the given by the following operator, 



Substituting Eq. fHSl) into the factorization condition and applying Eq.® we obtain, 

^2li?12 = 7^l2J■l2, (48) 

where the twisted i?-matrix, 7^i2, is the diagonal matrix, 

N N 

7^l2 = Xi ® X2 - 5^ ef ) ® e?) + Y ^ ® \ (49) 

i=l i=l 

Hence, applying unitarity the L = 2 factorization condition becomes, 

l<a2<ai<A'^l ^12 + Z^l<ai<a2<Af ^1 Xi (g) X2 

_ '7? ,<?N o("2a2) I Y^ 'P ^("lai) ^ „(a2a2) p 

- Z^l<ai<a2<A' '^1261 8)62 + Z^l<a2<ai<7V '^1261 ® 63 ii:i2. 

Using the following relation regarding the action of the 7^-matrix on elementary matrices. 



(50) 



„ (ii)^ U,) J «i(6,6)eS"^®e^^^^ for t = j 

-JZiiel ' ® e'^'-" = ■I , (51) 

I e\ ' 1^ 63 for i 7^ j 

we can immediately see that the ai ^ components of the summations in Eq.f l50l) cancel. Hence 
we are left with the expression, 

N N 

^el"") ®e("")i?i2 = 5^7^l2eS"") ® e^""\ (52) 

a=l 0=1 

which is true by inspection. 

We end by commenting that the form of J^u coincides exactly with the solution originally 
proposed in [15] for the rational SU{N) vertex models. This observation provides us a hint on 
how to proceed for arbitrary L. 

3.2 The general L case 

For general L we provide the following ansatz for the form of X'i...l, 

F1...L = M...L-^l...L, (53) 
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where the definition of Mi,,,l is given by the product of partial TV-matrices, 

M...L = Af2...LAfl,2...L = ^f{L-l)LJ^L-2,{L-l)L ■ ■ ■ Afl,2...L, (54) 

such that the partial A/'-matrices ^^(i+i)...^, are given by, 

A/'i,(i+i)...L = MiL-^fiiL-l) ■ --Uiii+i). (55) 

As with the L = 2 case, the terms J-'i,,,l coincide with the form of the solution given in [15]: 

j'i...L=Yl E (8)4""^^ll' (56) 

o-eSi, i<Oct(i)---"<t{l)<^ «=i 

where the symbol * in the sum f l56|) of ordered indices is to be over all non decreasing sequences of 
the indices aa-(i). The indices ^^-(i) satisfy one of the two inequalities for each pair of neighboring 
indices: 

a<T(i) < if cr(0 < o'(« + 1) 

(57) 

a<7{j) < if cr(«) > cr(« + !)• 



The specific choice on the second part of the sum ( 1561) and the form of the /2-matrix ensure 
that J^i,„L and Fi x are lower-triangular (see [T514T7] for details). Additionally, since each diagonal 
entry is non-zero, the inverse F^^^ is assured to exist. Note also that the tensor product term in 
Eq.( l56p is invariant if we apply the permutation to each index i, 

i=l i=l 

and hence an equivalent expression for the J-'-matrix is, 

^i-L=Yl E (58) 

We now are left with the task of verifying the factorization condition fl37j) . This involves a 
sequence of steps that we shall now detail. 
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3.2.1 Recasting the factorization condition 

In order to show the vahdity of the factorization condition (1371) for all a G S*/, we proceed 
much like [12] and take advantage of the decomposition of Fi i present in Eg. ( 15^ . In doing so 
we recast Eq. fl37|) as an equation involving J^i,,,l and the twisted i?-matrix. We then adapt the 
procedure first devised in [T5HT7] to tackle the problem. To this end we substitute Eq.f l53|) into 
Eq.(l371) to obtain, 

•^.(i..x)i?ii = Ar-} i)M...L-^i...L for all aeSL. (59) 

Recall that all elements of Sl possess a minimal decomposition in terms of adjacent permutations 
( 136|1 . For the expression A/'^J ^)A/i...l we offer the following result. 

Proposition 1. 

Kil...L)^i-L = T^\':!L for all aeSL, (60) 
where 'R-l'^^ follows the same decomposition rules as R\"\, 

where TZl'^^^P = 7^K^("i+i)}7^K2('^2+i)} _ _ _ (qi) 

and = Pa{a+l)T^a{a+l) , 

and lZi2 is explicitly given by 

Because the corresponding auxiliary operators 'R-'f^L provide a valid representation for Sl, 
the above result can be verified by showing that Eq. fl^U]) holds for only two permutations: the 
adjacent permutation ai2 and the cyclic permutation cxc = o'i2(y2z ■ ■ ■ For more details 

refer to Appendix C of [19j. 

Hence we recast the factorization condition in the following form, 

^.(i...L)i?.Jl = 7^^i^l...L for all aeSL. (62) 

We impose that the 7^-matrices follow the same left-handed and right-handed convention as the 
i?-matrices, 

'^l,2...,Ar = '7^1Af'^l(Af-l) • • • '^12, and 7ll...N-l,N = 'TIin712N ■ ■ ■'T^{N~1)N- (63) 
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Following the above convention the 7?.-matrices obey the same global unitarity condition as the 
i?-matrices, 

7^l,2...L7^2...L,l =Xi®---®Xi, (64) 
and being diagonal the 7?.-matrices commute amongst themselves. 

3.3 Verification of the factorization property 

We begin by applying Eq. (HO|) to Eq. (l62|) to obtain, 

R^i.:!^ = -^iTi.^tr (65) 

Since both R^^i^ and 7^1'^ ^ provide valid representations of Sl, we remark that Eq. (l65|) is in a form 
that one can readily decompose the permutation a. To illustrate this consider the permutation 
cr = (Jicr2 on the left-hand side and right-hand side of Eq. fl65|) respectively, 



^l...L - ^l...L ^l^.L _^ ^gg^ 

77-1 ^{(0-10-2) ^} 77 _ 77-1 'r,{^2 } 77 77-1 '73'f°"i ^ T 

where {ai, 0-2} e Sl- 

Since Sl can be constructed entirely from the adjacent permutations (Tj(j+i), j = (1, . . . , L — 1), 
we need only verify Eq.( l62i) for the adjacent permutation to guarantee its validity for all Sl- To 
this end we substitute a = Cj^+i) into the factorization condition (l62l) obtaining, 

•^f7,(,+i)(i...L)i?i(i+i) - ^i(i+i)-^i...L = 0, (67) 

where, 

apply Eq.l(4TJ 



77 D _ \^ \^ /0\^("<^(i)"<^(i)) d{'^} d'^jO + I) 

^<x,(,+i)(l...L)^j(j+l) - 2^ 2^ Vy%u+i)<^ii) ^a,(,+i){l...L)^l...L 

* L 

(tG5£ l<a<j{i)...Q!<j(L)<Af j=l 
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We now perform the change of variables aa{i) — )■ aa-^-_^_j^^a{i) to the summation indices Q;<j(i)- In 
doing so Eq. (!68l) becomes, 

** ^ I \ 

><T,(,+i)(i...L)ii:i(i+i) - 2^ 2^ Q9%(,+i)-w ^i-i 

where we have apphed the relabehng (Tj(j+i)cr = r of the elements of Sl- 

The summation EI<a,(i)...a,(i)<iv Eq.([69]) for is deceptively similar to 

^Y^x<a a. {L)<N Eq-dSH]) for in that its ordered indices are to be summed over all 

non decreasing sequences of the indices ar(i)- However the indices ar{i) satisfy one of the two 
inequalities for each pair of neighboring indices: 

ar{i) < ar{i+i) if a^{j+i)T{i) < a^{j+i)T{i + 1) ^^^^ 
ar(i) < ar{i+i) if crj{j+i)r{i) > (Tj(j+i)r(i + 1) 



Comparing fl57I) with fjTOj) . one can see that the only difference between them is the adjacent 
permutation crj{j+i) factor in the if conditions. For any given r G Sl, we focus on two integers, k 
and /, where T{k) = j and t{1) = j + 1. Using these two integers we examine how the elementary 
transposition 0'j{j+i) will affect the inequalities in (1701) compared to (!57|) . There are only two 
relevant cases, given by \k — l\ > 1 and \k — l\ = 1. When \k — l\ > 1 then the adjacent 
permutation crj(j_|_i) does not affect the sequence of indices a,- at all, meaning that ( |5711 and (ffOll 
are the same. To explicitly show that fl57j) and fITOl) are the same in this case we concentrate on 
the following four tuples, 

(«r(fc-l), «r(fc)) («r(fc), «t(A:+1)) («r(«-l) , «r(0 ) («r(0 ) C^rCi+l) ) , 

which are affected by the permutation crj(j+i). 

Focusing on the tuple (ar(fc-i), «t(A:)), there are two possible cases for the values of r(A; — 1) 
and T{k) given by, 

aj(j+i)r(/c - 1) > crj(j+i)r(A;) or (Tj(j+i)r(/c - 1) < o-j-(j+i)r(A;) (71) 
Since |fc — /| > 1 we remark that for both cases given in Eq. fITT]) : 
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• r(fc — 1) is invariant under the permutation cTj-q+i). 

• crj(j+i)r(A;) = <yj{j+i){j) = j + I - changing r(A;) to T{k) + 1 is not going to affect the values 
of the inequahties. 

An equivalently elementary analysis can be performed for the three remaining tuples. Hence with 
these observations, we are assured that the type of the inequality is unchanged by the addition of 
the permutation c^xi+i) " meaning that fl57j) and f lTOl) lead to the same results when \k ~ l\ > 1. 

We now concentrate on the case \k — l\ = 1. In this situation the permutation crj(j+i) does 
affect the inequality, leading to a difference between (!57|) and ( I70|) . To see this consider the case 
'^{^) = j ^-iid T{k + 1) = j + 1, where 

T{k) < T{k + 1) and 

^i(i+i)^(^) > + !)• 

Thus, what is usually a "<" inequality in f l57|) changes to a ">" inequality in (I70|) . Hence the 
summation in Eq. (l56|) . 

■■■'^T(k)'^T(k+i) aj<aj+i... 

changes when we consider the summation in Eq. fl69l) . 

E = E ■ 

■■■aT(fc)ar(fc+i) aj<aj+i... 

Equivalently, consider the case r(A;) = j + 1 and r(A; + 1) = j, where 

T{k) > T{k + 1) and 
(TjQ+i)T{k) < aj(^j+i)T{k + l). 
Thus, what is usually a ">" inequality in f l57|) changes to a "<" inequality in flTOj) . Hence the 
summation in Eq. (l56|) . 

•••aT(fc)"T(fc+i) aj+i<aj... 

changes when we consider the summation in Eq. (l69il . 



E = E ■ 

i^T(fc)"T(fc+i) aj+i<aj... 
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3.3.1 The case |A; - /| > 1 



Applying the above analysis we now consider the left-hand side of Eq. fl67p , 

and in particular, focus on the case where cr(fc) = j, cr(/) = j ' + 1 and |A; — /| > 1. In this case, 
^T<a^w-o^.iL)<N = ^I<a.(i)...a.(L)<W' ^^^^^e Eq.(^ becomes, 

(x,®x,^,-7^,(,^,)) EE E (73) 

|fc-i|>l CT(fc)=j,CT(i)=j + l 

Through inspection one can see that the condition |fc — /| > 1 means that aj ^ ctj+i in the 
summation - the easiest way to convince oneself of this statement is to consider the case \k — l\ =2 
(for clarity label / = k + 2) and look at the inequality conditions between aa-{k), C(a{k+i) and aa-{k+2)- 
Hence, through applying Eq. flSTl) . Eq. fl73l) becomes, 

apply Eq.JST} 



l<«^(l)...a^(i)<Af 



L 

{Q!cT(i)"<T(i)) d{o-} ^ g 
a{i) 1...L 

1=1 

3.3.2 The case \k - l\ = 1 

We now focus on the case (t(A;) = j, a{l) = j + 1 where \k — l\ = 1. In this case ^ 
Ka n\---a it\<n^ heuce we look at each term in Eq.(|72J) separately, 

L 

L 

1=1 

4 = 1 
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and, 



L 



L 

= '^3{3+l)T.i=lY.^^^^^-^^^lL^^^^^^YX<^^^ ^l!.L (76) 

1=1 

i=l 

We notice that the summations over the a^'s in Eqs. fTMl) and (!75|) are now of the same type as 
those in Eqs. (1761) and (177|) - i.e. they have only one "*" symbol. This is because the summation 
indices which cause the difference (a^ and aj+i) have been dealt with explicitly. 

We now subtract the aj ^ a^+i component of Eq. (!76l) from Eq. (!74l) to obtain, 

apply Eq. |(5T|| 

, ^ . 



<j('=)=j,<^(fc+i)=j+i 

L 



{Q!<T(0"<T(i)) D{fT} ^ g 
CT(j) 1...L 



1=1 



(78) 

Similarly we subtract Eq.(l77I) from the aj ^ Oj+i component of Eq.(l75l) to obtain. 



apply Eq. JSTIl 



<T(fc + l)=J,<T(fc)=j + l 



("£r(i)"<T(i)) p{(T} _ p. 



1=1 

i7^fc,fe + l 



(79) 

Finally, we now subtract the = Oj+i component of Eq. (1751) from the = aj+i component of 
Eq.( !75|) to obtain, 

L 



(80) 



-7? \^ S^* i<>l^("<^(i)"<^(i)) pW 

i=l 
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We proceed in Eq.( !80|) by first making the change r = aak{k+i) in permutation labels. Next we 
apply Eq. (HT]) to the permuted i?-matrix rI^\ to obtain, 

pM _ p{'^fc(fc+l)} p{o"} _ p tdW} /oiN 

^l...L - ^a{l...L) ^l...L - ^iO+l)^l...L- 

Taking into account the above change in permutation labels and decomposition of the permuted 
i?- matrix, Eq. llHOi) becomes, 

L 

i—l 

i^k,k+i ^g2^ 
^ V ' 

apply Eg.l lSTIl 

Hence we have verified the factorization condition for general L and A^. 



4 Twisted monodromy operators for = 3 



The purpose of this section is to show that the constructed F-matrices can effectively be used 
as similarity transformation in the simplest situation of the = 3 state vertex model. We shall 
present an algebraic derivation of the form of some relevant elements of the monodromy matrix 
in the F-basis for arbitrary weights. Here we represent the = 3 monodromy matrix as, 

( Af^i,) aS(^) B^lM \ 

ra,i...M= . (83) 

V ci'!M c['!m d^...m I ^ 

In general, a given monodromy matrix element Xi...l(/^) can be transformed to a new operator 
-^i...l(/^) having a much simpler quasilocal form with the help of the F-matrices [1]. Since Fi i 
and J-'i...L are related via the multiplication of a diagonal matrix (l53!l . it is enough to compute the 
following non-trivial twisted operators. 



Xi...l(/x) = j'i...LXi„,i(/i)j'i:.!i. 



M) 
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In what follows we shall present a conjecture for the expressions of the twisted monodromy 
operators Di...l{ij)-, C[^\(ijl), B^\{ijl), C[^\{ijl) and B^\(ijl). Our results are built up from an 
analysis performed in the cases L = 2 and L = 3 relying only on the identities derived from 
the Yang-Baxter and unitarity relations. Fortunately, this study is sufficient to foresee the main 
structure of the mentioned twisted operators for arbitrary L and general Boltzmann weights 
without relying on any specific parameterizations. 

We now list the final expressions for the above mentioned twisted operators, 

L 

Di...M = diag {631 (/^, ^i), bsiil^, ^i), asil^, ^i)}^ (85) 

i=l 

C?IM = X]c32(//, diag |62i(/^, ^i), ' "^(/^^ CO^afe, 6)| (86) 

BflM - Y.c,,{n, i,)ef^ (g) diag ^,),h2{^, ^i)e2{^i, ^ b3^(£'^^^g3f^^, } (87) 



L 



, C3I (/^, 61 )^32(^, 62)^21 (61 , 62) ^(23) „ ^(12) 

^.4--. M6::^;j (^8) 



il5^'2 



^ / b2l{l^,^i) &32(/^, 6)^2(6, 62) / CA/J/t c ^\ 

, Ci3(/i,{iJ632(/i,62)ci2(6i,62) (32) „ (21) 

«l#i2 

f, / tN/J^t t\ ^32(/i,Ci)6'2(C/i,Ci) a3(/^,^j) 



i=l 



?^2i(ei, 62)^2(62, Ci) ' hiiii, 61)^3(61, J , 
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where the auxihary functions Oi{^j,^k) are given by, 

{O'ii^i^ ^k) for j < k 
z = 1,2,3. (90) 
1 for j > k 

We begin by exphcitly verifying that the L = 2 case holds for each of the twisted operators 
given below. 

4.1 The L=2 case 

In this section we shall calculate the twisted operators -Di2(/i), C^{fi), B^{fi), C[2\^J') and 
(h) directly from the similarity transform ( 18^ . We shall then detail the necessary Yang-Baxter 
and unitary relations for the entries obtained from Eg. (18^ to match the corresponding entries in 
the L = 2 operators given by Eqs. (18511891) . The expressions derived directly from Eq. (!84l) have the 
following structure, 

DuiiJ,) = diag {63i(/i, ^i), &32(/i, 6), asl/i, 6)}i ® diag {631 (/i, 6), &32(/i, 6), asl/^, 6)}2 



(91) 



(92) 



All the K entries of the above matrices do not immediately match the corresponding entries 
calculated from the L = 2 expressions of Eqs. fl85|) -f l89|) . It is possible to simplify these entries 
using only certain relations coming from the Yang-Baxter ([T]) and unitarity (j2]) relations. The 
technical details are quite cumbersome and thus have been deferred to Appendix B. In Table ([T]) 
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Operator 


Y-B equations 


Unitarity equations 




(II2])-{3,1,2} ([MD-{3,l}-{3,2} 






(II1-{2,3,1} (Il6])-{3,2} (Il9])-{3,2} (I2l])-{3,2,1} 


(IIOD-{3,2} 


^(2) 
^12 


(II1-{3,2,1} (inD-{2,3} (Il8])-{2,3} (I2l])-{2,3,1} 


(IIOD-{2,3} 


'-"12 


([I3])-{1,3,2} ([16])-{3,1} ([I9D-{3,1} (!20])-{3,2,l} 


(!9i)-{2,i} m-w) 


^12 


([I3])-{1,3,2} ([I7D-{1,3} ([T8D-{1,3} (!2I])-{1,3,2} 


(!91)-{1,2} ([I0D-{2,1}-{3,1} 



Table 1: Required Yang-Baxter and unitarity relations for L = 2. 



we provide a summary of the unitarity and Yang-Baxter relations that are required to simplify 
such non-trivial entries in order to bring the operators in the form given by Eqs. fl85p - (|5^ . 

We remark that in this table a given algebraic relation among weights is referred to the equation 
number together with its respective indices {i}, or fc}. For example, the symbol ([9])- 

{1,3} refers to the equation ^13(^1, 6)^3i(6, ^i) + £13(^1, 6)c3i(6, ^1) = 1 whilst ([T2D-{1, 3, 2} 
means the relation 613(^1, ■^2)^12(^5 'C3) = ^i2('Ci) ^2)^i3('Ci) ^3)- We observe that there exists obvious 
equivalences for some equations such as for Eq. (fT2|) the indices A;} or {i,/c,j} leads to the 
same relation. 

At this point we remark that the L = 2 case does not allow one to verify the entries of the 
diagonal matrices of the double summation from C*j^^^ and b[^\, we have only verified the coeffi- 
cients of the double summation terms, given by ^3i(M.g^i)'>32(M.g^.c.ite,,g.,) c,3(/.,g,,)te(A..gyci.te„g,,) 

respectively. In order to obtain the entries of the aforementioned diagonal matrices we need to 
consider the L = 3 case. 



4.2 The case L=3 

We now consider explicitly calculating the L = 3 case using Eq. (l8^ . Here it is only necessary 
to consider the operators c[^^{^) and -8(23 (/^) compute the entries that are completely missed 
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by the L = 2 case. We begin with the operator CjgslAt) whose entries are proportional to, 

As with the L = 2 case, we shall detail the Yang-Baxter and unitary relations necessary for the 
entries obtained from Eg. f l84p to match the L = 3 case of Eg. (188 p . In Table we provide a 
summary of the unitarity and Yang-Baxter relations that are reguired to simplify all the entries 
in order to bring the operators in the form given by Eg.f l88p . The technicalities of this calculation 
are rather involved and thus have been deferred to Appendix C. 

We now turn our attention to the operator 5(23 (/i) whose entries are, completely missed by 
the L = 2 case are given by, 

ef) ® ef ) d^ag{a[''^\ «f ^^js + ef^ ® d^ag{(3[''^\ /3f (3^''% ® ef ^ 
+ d^ag{^r, 7^ 7^ ® ef^ ef^ + ef^ ef ^ d^ag{6r, 4'% 
+ ef'^ ® dzag{<f>['''\ <pf^\ <pf'\ ® ef ^ + dtag{uj[''^\J^-\ ujf^\ ® ef ^ ® ef\ 

The technical details entering in the simplifications of these entries are similar to that conducted 
for the operator 6*^23 (/i)- We therefore restrict ourselves in presenting only the reguired Yang- 
Baxter and unitarity relations that are necessary for the raw entries obtained from the L = 3 case 
of Eg.f l84p to become the corresponding L = 3 entries in Eg. (!89|) . This is summarized in Table 

4.3 The general L case 

At present we just have a proof for arbitrary L in the case of simplest the twisted operator 
We now present the proof of Eg. fl85p for general L. Our verification is an adaptation 
of an argument first given in [15] for the eguivalent operator. To begin we note that the operator 
Di...l{^i) can be expressed through, 

ef'^D,,,,M = ef')%,,...MeT\ (94) 
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Entry 


Y-B equations 


Unitarity equations 


(Ci) 

i 


([H-{3,2,1} (I20])-|3,2,1} 


(Iini)-{3,2} 




(1201)-|3,2,1} 


(Iini)-|3,2} 


o 


(EOD- {3,2,1} 


(iini)-{3,2} 


1 i 


([H-{3,2,1} (I20])-{3,2,1} 


(IT0])-{3,2} 




([Il-{3,2} (I20])-{3,2,1} 

\! y L ' J v y L ■ ■ J 


(iini)-{3,2} 




(II1-|3,2} ([l5])-{3,2} ([l6])-{3,2} (I20])-{3,2,1} 


(Iini)-|3,2} 


/ 1 


([I3])-{3,2,1} (II9])-{2,1}-{3,1} (I20])-|3,2,1} 


([IO])-{2,l}-{3,l}-{3,2} 


(Ci) 


([I3D-{1,3,2} (II6])-{2,1} (Il9])-{3,2} (I20])-{3,2,1} 


(iini)-{3,2} 


7f^ 


([I2D-|3,2,1| ([I5D-|3,2| (II6D-|3,1| (|20l)-|3.2.1} 


(fT0l)-|3.2} 

y L * J 


I 


([I3D-|3,2,1| 




z 














m-{w] m-w^] m-w) 


(^-{2,1} 
















m-w^] m-{w] (ii9])-{2,i}-{3,i} m-i'^M 


([I0])-{2,1}-{3,1} 




([T3l)-{1,3,2} ([I6D-{2,1} (I20])-{3,2,1} 


([ini)-{3,2} 




([I2D-{3,2,1} ([Il-{3,1} m-{W] 





Table 2: Required Yang-Baxter and unitarity relations for the entries of C 
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Entry 


Y-B equations 


Unitarity equations 


al ' 


(II3D-{2,3,1} ([I8])-{1,2}-{1,3} (I2I])-{1,3,2} 


©-{1} 




ffT3l)-|l 3 21 ffT71)-il 21 (lT8l)-i2 31 fl2Tl)-il 3 2| 


(I8l)-i2| ffT0D-|2 11 


(Bi) 


ffT2l) - 13 1 2|(|T3l) - 12 3 11 
(HD - |3,2|([I7D - il,3|(l2lD - |1,3,2| 


(IIU])-{2,l}-{3,l}-{3,2} 




([13D-|2,3,1| fl2T])-i 1.3.21 


f[T0D-i3.2l 




([18l)-i2,3j fl2TD-il.3.2l 






f[T2D - |3, 1, 2|([T3D - {2, 3, 1 j([Tl - |3, 2| 
ffT5D - |3,2|([I7D - |2,3|(I2I1) - |1,3,2| 


([IQD-{2,1} 


(Bi) 
/ 1 


(12I1)-|1,3,2| 




(fii) 

/ 2 


(I2I1)-|1,3,2| 




/3 


^-{3,1,2} (EB-1 1,3,2} 




i 


([I2|)-|3,l,2} ([I81)-|1,2}-|1,3} (EB-1 1,3,2} 


dHD-ll} (fT0D-i2.1j 




([I3l)-|l,3,2} ([I4|)-|2,l} (II71)-|1,2} (I2l])-|l,3,2} 

yi j' L * * J yi j' L ' J \i y L * J yj y L ' ' J 


(IH1)-|2} (fT0D-i2.1| 




(II7D-{1,3} 

y! y L ' J 


(fT0D-|2.1}-{3.1} 




(II2D-{3,1,2} 






(II1-{2,1} 


(iinD-{2,i} 






(iinD-{2,i} 


^1 




([IQD-{2,1} 






([IQD-{2,1} 




([I2D-{3,1,2} 


([I0D-{2,1} 



Table 3: Required Yang-Baxter and unitarity relations for the entries of B 
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where the Weyl matrices ei^^^ project out the operator Di from the 3x3 matrix expression 

in auxihary space Aa- 

Using the above expression we now consider the action of the J^-matrix on 



apply Ea.BSt 



,(33) X- n {,A - sr sr* /c?\„(°-w"-w)„(33) pW -r (.aJ^^) 

3a J-1...lI^1...L[N - Z^aeSi Z^l<a,n)...".(i,)<3^^^<7(i) ^l...L ' a,l...L[N ^'a ' 



i=i 95 

L 

— Z^ae5L ^l<aa(i)-a<T(L)<3 VCy^fT(i) ' a,a{l...L){H')fia ^l...L- 

i=l 

In what follows, we separate the sum over the indices acr{i) according to the number of occurrences 
where ao-(i) = 3, 

L *{k) L 



(t65l fc=0 l<aCT(i)...aCT(L)<3 «=! 

where, 



*(fe) * 

E = E • (97) 

l<aCT(l)---a<T(L)<3 l<a<T(l)---<»CT(L)<3 

(a<T(l)'-"'"<T(L-fc))e{l'2} = ("CT(L-fc + l)'"""<T(L))=3 

L 

Hence the tensor product of the Weyl matrices (^e^""^''-'"'^*'-'^ for each value of the index k becomes, 

i=l 

L L-k L 

^(ao^(i)a<T(i)) ^ ^^(atr(OQ<T(i)) ^ ^(33) 

i=l i=l i=L—k+l 

leading to the following form for Eq. fl^B]) . 

L 

(33) T- n (,,\ — sr^ Y^*(fc) /<^(^<^{i)'^'y{i)) 



Ga J^l...LDi,„L{fJ') — Zl<je5i Z]fc=0 Z]l<a^(i)...a^(i)<3 Qy'^a(i) 

i=l 

(33) (33) n / \ (33) T-> ( \'-r ( \ (33) 

XCa e^(^)i?aa(L)(/i) . . . e^(^_;,^^)i?aa(L-fe+l) (^)7a,a(l...(L-fe)) (/^jCa i?!,,,^. 

(98) 

In the above expression we realize that the action of the Weyl operators on the i?-matrices 
leads to, 

L L 

j=L-fc+l j=L-fc+l 
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and hence Eq. (p8|) becomes, 



xe 



L 

("<T(L-fc)Q<T(L-fc)) o ^(q<,(i)Q,(i)) p {,,\J33) ^ J33) 

a{L-k) 



j=L-k+l 

Since ao-(i) G {1,2}, i = {1, . . . ,L — k), in the above expression we reahze that the action of the 
Weyl operators on the /^-matrices leads to the following simplified expression. 

The above considerations lead the following result. 



L \r^*{k) y-rL—k 
k=0 ^l<o<^(i)...aCT(L)<3 •!■ •l-j= 

nt=L-.+i«3(/i,U))(8)e^"'^i?l: 

diag {631 (/i, ^j), 632(/i, 6), a3(/^, 6)}^ -^i.-.l, 



J^l...LDi,„L{fi) - Eae5z, Zlfc=oI^l<«,(i)...a<,(i)<3ni=l ^3a,(,) (/i, ^f7(j) ) 

L 



ultimately verifying Eq. flHSl) for general L. 

We do not have a proof for general L in the cases of the twisted B and C operators. The 
necessary underlying recurrence relations to carry out such demonstrations have thus far eluded 
us. However, with the help of the Yang-Baxter solution presented in Section [2] we have been able 
to verify Eqs. fl86ti89p for L = 4 explicitly. Considering the generality of such solution, this is strong 
evidence supporting our conjectured expressions Eqs.( !86ll89l) for general L. 



5 Basic domain wall partition functions 

The purpose of this section is to start the formulation to compute certain domain wall partition 
functions (DWPF's) associated to the = 3 vertex model with arbitrary Boltzmann weights. 
The domain wall boundary conditions correspond to certain fixed statistical configurations for the 
horizontal and vertical edges at the top and bottom of the square lattice, see for instance [911261 127]. 
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From an algebraic perspective these objects can be expressed in terms of the expectation values 
of combinations of the creation and annihilation operators B[^\{jj) and C[^\{h') on some pseudo- 
vacuum states P, El EH [29]. In Figures [2] and |3] we have depicted the graphical representation of 



such domain wall partition functions. 



1^1 



I'M 



M 



Figure 2: Graphical representation of a mixed DWPF of type C. 

In order to compute such partition functions we need to first define the corresponding reference 
states as well as their properties under the action of the F-matrix. The = 3 state vertex model 
has three possible pseudovacuum states which we label as follows, 

fi\ fo\ fo\ 







|2). 



1 



|3). 







(99) 



for a = 1, . . . ,L while the corresponding transpose states are denoted by, 



,(1| = (1,0,0), , «(2| = (0,1,0), , ,(3| = (0,0,1),. 



(100) 



We also use the following convenient notational conventions with the outer products of reference 
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M ■ ■ ■ ^1 

Figure 3: Graphical representation of a mixed DWPF of type B. 

states, 

p 

\i)ar...a^ = (^\i)a, i = (1,2,3), 

i=i 

and similarly with the transpose reference states. 

One then finds that the action of a Weyl basis element on a reference state and a transpose 
reference state respectively is given by. 



Using the above observation we give the following results: 



a{k\Sij 



(101) 



Proposition 2. 



J^l...L\i)l...L = \i)l...L , l...L{i\J^l...L = l...L{i\ 
J^llL\i)l-L = \i)l...L , l...L{i\J^i,]i 



(102) 



1...LW 



Proof. From Eq.f llOip we notice that the only element of the expression (3^6^"^"^'' in Eq.( l53|l 



which acts non trivially on the reference state \i)i,,,L and its transpose is when ai = 0^2 = 

■ ■ ■ = ai = i. From our detailed analysis of Eg. (155]) in Section we know that this particular 
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set of a values in the sum of Eq.( l57p only occurs when o" = X is the trivial permutation. Since 
I^Y^i^ = Xi (g) ■ ■ • ®Xx, this proves the first two relations in Eq. fll02p . 

Focusing on the final two results, from the above analysis we find that the 1st, ^^^y^th and 
3^th row and column of T\,,,l are given explicitly by 1 in the respective diagonal position, and 
zero everywhere else. Noting the following cofactor values concerning the 1st row and column of 

-^l...L, 

det {(J'i...L)fci}fc=2...3i = ^Qi{Tx...L} 



1=2. ..3 



L 



det{(J'i...L)fcJfe=i,.,(,_i)(,+i)...3L = for j G {2,...,3^}, 



det{(J-i...L),J ,=2...3i = forzG {2,...,3^} 

i=l...(i-l)(i+l)...3-£' 

)...3i — U iUl J C . . . , 

1=2. 

and similarly equivalent values for the ^^^th and 3^th row and column of 
We now recall Cramer's rule for the inverse of a matrix, 



(-^i-'^)^. = Hpf/r\ ^^^ {(^i...l),J.^....-i.+i...3. . (103) 
■> aex \j- 1,,, Li i=i...i-i,i+i...3i 

Hence applying this equation we obtain, 

(^i:.'L)n = 1 

(^iT.y.i = forzG{2,...,3^} 

{^i'l),, = forjG{2,...,3^}, 

and similarly equivalent results for the ^^^th and 3^th row and column of J^i,„l - thus verifying 
the last two results of this proposition. □ 

We now start by computing the four possible basic non-trivial DWPF's consisting of the 
expectation values of tensor product of operators carrying a single charge index. These are the 
building blocks required to construct the domain wall partition functions involving mixed fields 
such as Ci^.?i(z/) and Cf^Liu) or B^^ M and ^^^(/i). 

5.1 Single DWPF for C^^) 

We now consider our first basic DWPF, labeled Z},"' ^ which is given explicitly by, 

Zf^'\{y],m = i...L(2|Cf)jz.^)...Cf)jz.i)|3)i..x. (104) 
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Applying the J-'-matrix similarity transform to each of the C^^-* operators in the above expression, 
and using the results of Eq. fll02p . Eg. (11041) immediately becomes, 

Zf''\{y].m = i...L(2|C'i^.^^K)...Cf.)^(z.i)|3)i..x. (105) 

We now are going to obtain a recurrence relation for the above expression. This is done by 
inserting a complete set of states between the operators C\ ^{1^2) and ^.(^i) ^° obtain, 

Zr''\{^},m) = E^=ii...L(2|C'f.)j^^x)...C'f.)jz/2)|2),|3)i . 

x,(2|. .(3|Cf.)jz/0|3)i...L. 

It is elementary to show that only the above terms of the complete set of states are non zero by 
considering the action of e^^^^ on the reference states. Through the property p(2| i...L(3|ep^^ = 
i...L(3|5pi, one is able to derive the following identity. 



L 

.(2) 



p(2| ^...l{3\C^^\{u,)\3)i...l = cs2{lyuQUas{uu^^)es{^^,Q. (107) 

We now focus on the expression i j^{2\C[^\{i'l) . . . C[^\{u2)\2)p\3) i...l . Using the fact that 
ep^'^''|2)p = 0, we can discard the value p in the summations of the C[^\{h') operators. Hence 
applying this fact, and elementary matrix multiplication, the expression now becomes. 



,...i(2|C|'>J..i)..,C|'>Jt.,)|2),|3>.^^. 



(lOJ 



where, 

L 



C'^fliy) = c,2{v,ii)ef^^diag (&2i(^,e.), , il'f.tl , v <^^(^^^^)H^^.^i) \ ■ (109) 

~f I 032[Q,C,i)ty2{i^i,Q) )i 

i^l,p 

We now offer the following comments on Eq.f llOSp . Firstly, the product 11^=2 fe32(g!'^g''?ef(g g;.) 
is independent of the value of I2, . . . Jl due to the condition h ^ ■ ■ ■ II ^ V- Secondly, the 

operators C\,],i {v) act trivially on the vector space V^, meaning that we can decrease the number of 

'^p 
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relevant vector spaces in the reference states by one. Thirdly, when the aforementioned coefficient 
is taken out of the sum over I2, . . . the following simplification occurs, 

1...^ {2\C^':'^\ul) . . . C?.:'l\u2)\3) ..... = ..... (:i\&?„, {vl) • • • C'l'L (z^2)|3) 1.... 
Hence Eq.f llOGp becomes the following recurrence relation, 

zrm,m = E,i,c3.(.„e,)ni;^i^nt^.&3.(-„e.) ^^^^^ 

5.1.1 Exact solution 

We now can apply induction to verify the complete algebraic expression for the C*^^-* type 
DWPF in terms of Boltzmann weights. To begin for L = 1 and L = 2 we have, 

Zf''^\^l,^l) = C32(z/i,^i) 

7(C>2)/r T rp1\ _ C32 (i/l ,gl )C32 {U2 ,6 )a3 {1^1,^2)93 (6 ,gl )fe32 (i^2 ,6 ) 

I C32 {ui ,g2 )C32 (i^l ,g2 )a3 (t^l ,gl )^3 (gl ,^2)^32 (i^2 ,^2) 
632(51,6)62(6,51) 

We now offer the following general result. 
Proposition 3. 

z^'^'^VX;/! — rr^ r (1,^,^71 °3(i^j,ga(fc))b32(i^fc,5a(j))e3(ga(fc),5,T(j)) (^^^) 

Zl - 2.ag5^1L=lC32(^^^,4<x{^)jlll<J<fc<L 632(C.w,?.(,))02(5.(,),?.(fc)) ^^^^^ 

Proof. Noting that the above formula is correct for L = 1, 2, we assume that it holds for some 
L, and focus on the L + 1 case of Eq.f lllOp . 



y{C,2) _ sr^L+l / N rri+l a3(t^l,$m)g3(gm,gp) Tf U ( 



L+1 „ f N TT "3 {l^j ,^a(k))h2 jl^k ,icT(j) )^3 (6(fc) ) 



X E^gs^f) ni=2 C32{j^t,^a{t)) n2<i<fc<L+l 



''32 (6{fc) )^2 (?<T(j) ,?CT(fc)) 



5 



where the sum over the permutations with superscript p is given by, 

L+l L+1 L+1 

^ = ^^--- ^ for ^ (T3 ^ ■■■ ^ a^+i. 

'^^^L ^P itp 

The verification of the proposition follows immediately through the change in label, j9 — 5- (Xi. □ 
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5.2 Single DWPF for 5^^) 

We now consider the equivalent DWPF expression for the 5) operators which is given 
exphcitly by, 

= ^...Lm^?..M...Bf}M\'^)i...L ^^^2) 
= i...L(3|5;^.)j/ii)...5f)^Mi2)i...,., 

where we have apphed the results of Eq.( ]102p to twist the operators. As with the previous DWPF, 
we now obtain a recurrence relation for the above expression. We insert a complete set of states 

~ (2^ ~ ('2) 

between the operators B\ and B\ j^{^2) to obtain, 

^,{2\. .{?,\Bf}M---B?^..M\2),...L. 

Once again one can show that only the above terms of the complete set of states are non 
zero by considering the action of e*^^^-* on the reference states. Taking into account the property 
ep^"* |2)p|3) = 5pi|3)i...i one finds that. 



~ (2) ~ (2) 

Let us now consider the expression p(2| i...l {3\B\ i^{n-2) ■ ■ ■ B\ ^{^,l)\2)i,,,l. Using the fact that 
we can discard the value p in the summations of the B^\{jj,) operators. Applying 
this fact, and elementary matrix multiplication, the expression becomes, 

p(2|. .(3|5f)j/i2)...5f.)j/XL)|2)i...^ 

^ ^ (114) 

= E'^'^.....^! nt=2&32(/i„ep)^2to„ep)p(2| . .(3|5S)(/i2)...5S^(/iL)|2)l...L, 

where, 

L 



I ^32(6, 6)^3(6, 6 



As before the product 11^=2 ^32(/^i) ^p)^2(^«j , ^p) is independent of the value of ^2, • • • , II and the 
operators B\'\^ (fi) act trivially on the vector space Vp, meaning that we can decrease the number 
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of relevant vector spaces in the reference states by one. Applying these simplifications we find 
that, 

mf:'l\fi2) . . . Bt'i\fiL)\2) 1.... = 1.... {3\BfL . . . BfL (/^l)|2) 1.... 

i2/'-'7^'L/P 5^P ^P /P #P 

Consequently, Eq. f lllSp becomes the following recurrence relation. 



x4^jf({/i},{oi/ii,g- 

5.2.1 Exact solution 

We now offer the following general result, 

where we note that the verification of Eq. (11160 follows exactly from the verification of Eq. f illip . 

In what follows we shall consider the basic DWPF's that are constructed from considerably 
more complicated twisted operators. Nevertheless, the task of obtaining the explicit forms of the 
DWPF's of type and B^^^ poses no greater challenge than what we have experienced thus far. 

5.3 Single DWPF for C^^) 

The DWPF expression for the C'^\ operators is given explicitly by, 

zf^'\{y]Ai}) = i...L(i|C^i!.^z.K)...C7i^.^J^i)|3)i...L ^^^^^ 

= i...l(1|C'«Jz.^)...C';^)Jz.0|3)i...,, 

where we have applied the results of Eq. (1102^ to twist the operators. 

Before inserting a complete set of states between operators we shall first consider the expression, 
^i^^l(^i)|3)i...L) s-iid in particular we notice that the action of the matrices e'^^'' ® e['^^^ on the 
reference states |3)/^ 13)^2 is zero. Hence if we insert a complete set of states in between the 



(115) 
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operators C[^\{i'2) and Cl^^lui) we obtain, 

zr''\{^},m) = e^=ii...l(iic';^.^j^l)...^;^.^j^2)ii)p|3). . 

(118) 

x,(l|. .(3|C;^.^^(z/i)|3)i...^, 

which is very similar to the corresponding expression for C^'^\ 

We now use the identity p(l| {3\e\ = i...L(3|5pz and by elementary matrix multiplication 
we obtain, 

L 

i....(3|c7i!.^jz/i)|3)i...L = c3i(i/i,ep)n«3K,6)^3te,ep)- 



i = l 

7^P 



Focusing on the expression i,,,l{1\C[^\{i'l) ■ ■ ■ C[^\{i'2)\'^)p\^) ^---l we note that due to the 
elementary relations in Eg. (110 II) . there is no possibility of the Wely matrices ej^ ® ej^ in the 
expressions for the C[^\ operators to produce anything but zero, hence they can be discarded 
from the calculations. Additionally, using the fact that e^p^^\l)p = 0, we can discard the value p 
in the summations of the C*|^^^(z/) operators. Applying the above facts, and elementary matrix 
multiplication, the expression now becomes. 



i...L(l|C'«,(^^,)...C'i^),(z/2)|l)p|3)i . 



:ii9) 



where, 

C'.':.'i(^) -C3i(..;,)e,'"'(g)dfaj ''f 'ff . (120) 

~i lc'2i(4/,4i)c'i(4i,4/) 032(4«,4j) Ji 

We note that the product 11^=2 621 (g^^^g^Tef^g g,.) independent of the indices 12,...,Il, and 
the operators C^i,''^l (ly) act trivially on the vector space Vp, allowing us to decrease the number of 

't^p 

relevant vector spaces in the reference states by one. Applying such simplifications one finds the 
relation, 

E\......^i 1.- {i\c['i'\i^L) . . . c^l:'l\i^2m ...L = .... {i\c['l (z/^) . . . c^, (z/2)|3) .... 
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Again, we can discard the terms involving ef^'' ® e[^'^^ in the above C^^},l operators because they 



always produce zero. As a consequence of that, Eq.f lllSp becomes the following recurrence relation. 

The same type of arguments presented in Section 15.1.11 can be used to provide the general 
solution to the above recurrence relation. It is given explicitly by, 

5.4 Single DWPF for type B'^^^ 

We now consider the equivalent DWPF expression for the operators which given explicitly 

by, 

Z^L''\{f^},{0) = i...L(3|i?;^,)j/ii)...5^)jML)|l)i...L ^^23) 
= i...L(3|S;^)j/ii)...5«J^i)|l)i...z., 
where we have applied the results of Eq. fll02p to twist the operators. 

As with the case for C^^\ before we insert a complete set of states between operators we shall 
first consider the expression i i{3\B^\{fii), and in particular we notice that the action of the 
matrices e^"^^^ ® e;^^"* on the transpose reference states ij(3| (g) /2(3| is zero. Hence if we insert a 
complete set of states in between the operators B^li^{fii) and B^li^{n2) we obtain, 

^''''^(M,^) = E^=ii...L(3i5;^.^j/ii)ii)pi3)-. ^ ^ 

7^p (124) 
x,(1|..^..(3|5«^(/X2)...5;^.!^(/xl)|1)i...l, 



Using the fact that e] |l)p|3)i...L = Spi\3)i,„L one is bale to write the following identity, 

i...z.(3|EWj/iO|l)p|3)^^.. =c,3(/ii,ep)n , (.""f^nf, ,Y 

In parallel to Section [531 due to the elementary relations in Eq.f llOip . there is no possibility 

(21) (32) ~ (1) 

of the Wely matrices e^^ e^^ in the expression for the B^ ^ operators to produce anything 
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but zero, hence they can be discarded from the calculations. In addition, using the fact that 
Ijcp^^'' = 0, we can discard the value p in the summations of the Bf^j^ operators. Considering 



p 



these facts together with elementary matrix multiplication, the expression now becomes, 



.(3|5l^.^^(/i2)...i?}:.:j/iL)|l)l...L 



where. 



(125) 



Note that the product 11^=2 ^3i(/^i' ^p)^i(^'j' ^p) (^^es not depend on the indices /2, . . . ,/l and 



the operators -Bi^.'.l act trivially on the vector space Vp, allowing us to decrease the number of 



relevant vector spaces in the reference states by one. These simplifications lead to, 

1.- m^'±\fi2) ■ ■ ■ B^':'l:\fiL)\i) 1.- = 1..- (3|5«, (/X2) . . . B^'L (/iL)|i) i..- 

'27^'-'7^'l7^P ^J' J^P itp J^P 

f21) (32) ~ fl) 

Again, we can discard the terms involving e^^ ® in the above -Bi. .l operators because they 
always produce zero. Hence Eg. (11241) becomes the following recurrence relation, 

4"^^(M,{0) = E,iic,3(/.,,e.)ni^fiSiSnt2&3i(/..,e^ ^^^^^ 

x4^-^({/i}, {01^1,4)- 

The exact solution to the above recurrence relation is, 

iif^fd^}) -l^.esAU=lCl^{f^i,U^))nl<j<k<L fe3i«.(fe),C.o))e3{C.o),«.(fc,) • 



6 Mixed DWPF 

We start by discussing the preliminary steps to compute the DWPF's involving mixed opera- 
tors. The elements of the Bethe state vectors built from C fields for the fifteen-vertex model are 
given by, 

Cl^lMCtk^,.,) . . .cS(z.0|3)i...z., (129) 
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where ii, . . . ,ii G {1, 2}. Here we are considering the vector |3)i..x as our starting ferromagnetic 
reference state. 

We now introduce the integer M, M < L, which indicates how many type "(1)" operators we 
have in our state vector element. Given M, every element of Eg. (11291) will generally consist of an 
expression of C^^^ and C*^^^ operators in no particular order. Applying the following Yang-Baxter 
algebra expression generated from Eq. (jlj) , 

021 021 

it is possible to commute all the C^^^ operators to the left, leading to the general expression, 

ctl{i^L)ctl{u,^,) . ..ctki^.m,..., ^^^^^ 

where the coefficient 4>^M'{a}^^ ^ which is constructed from the Boltzmann weights 03, 621 and C12, 
generally depends on the initial value of the indices ^l, . . . ,iL, the integer M and the particular 
permutation a. 

From Eq. (I13ip we have that the most fundamental mixed DWPF expression of type C will be 
of the form, 

(132) 

where for clarity we assume that qi < q2 < ■ ■ ■ < Qm ■ We briefly note that the configuration in 
Figure [2] corresponds to Eq. (11321) for qi = 1, q2 = 2, . . . , qm = M. 

Similarly for B, the elements of the transpose Bethe state vectors are given by, 

i...^(3|i?^l(/.05!l(/i2) . ..Btlif^L), (133) 

where again ii, . . . ,iL G {1, 2}. 

Using M, M < L, to indicate how many "(1)" operators we have in our transpose state vector 
element, we apply the following Yang-Baxter algebra expression generated from Eq.(jl]), 

B'tlMB'^l.iu) = ^(/., u)B^l,{u)B^;'lM - ^(/., u)b['ImB^;'1,{u), (134) 
O21 O21 
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to commute all the B^'^^ operators to the right, leading to the expression, 

= EaGSi ^Af,{ar^ l...L(3|S[^.^i(/ia(l)) • • • ^L.X (/^-(M) ) ^^.^ ) • • • Sj^.^^ (/ia(L) ) , 

where the coefficient i^^^ia}^^^ constructed from the Boltzmann weights 03, 621 and C21 and 
parallels 4>M{a-}^'' from Eq. fll3ip in structure. 



From Eq. fllSSp we see that the most fundamental mixed DWPF expression of type B will be 
of the form, 

#91-. -gM 

(136) 

For clarity we briefly note that the configuration in Figure |3] corresponds to Eq. ( I136P for qi = 
l,q2 = 2,...,qM = M. 

We now devote separate subsections for the explicit evaluation of Eqs. f ll32p and f ll36p . 

6.1 Mixed DWPF for B 

Applying the J-"- matrix to Eq. fll36p to twist the monodromy operators we obtain, 

45/(M, m = i...l{^\b['1m . ..B'^^AmWilLim+i) ■ ..Bf^Af^L):F,...L\i),,...,,^2) ..^ . 

(137) 

We note that there is an J^i,,,l operator on the far right of Eq. (11370 as there is no equivalent 
expression such as Eq. (11021) for mixed reference states. We proceed by inserting two complete sets 
of states, but for clarity we shall do this in separate stages. Consider ffist placing a complete set 
of states in between the operators B\^\{fiM) and B^\{fiM+i) to obtain, 

zi%m,m = Ei<p,<...<p„<L i...L(3|5«j/ii)...5;^.^j^M)ii),,...p„i3) .... 



X pi.-.pM 



;i| 1.... (3|5f.).J^M+l)...5f.^Jy^L)-^l...L|l)......M|2) 1.- . 



(138) 

In order to verify that the above terms of the complete set of states are the only ones which are 
non zero, consider the expression, 

l...L(3|5«.J/il)...5«.J/iM), 
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and recall from Section that only the Weyl matrix ef^ (as opposed to ef^ e^^) gives non zero 
terms when applied to the reference state i...l(3|. Hence since i(3|ep = this means the only 
terms in the complete set of states which are non zero are, 



2_. |1)pi-pm|3)_^ i-i pi...pm(1L i-'f' (3|- 

l<Pl<-<pjV/<i' 



We now insert another complete set of states in between B^\{jj,L) and J^i,,,l in Eq. fll38p to obtain, 

zSLm,m) = Ei<p,<...<p,,<L i...L(3|5«.j/ii)...5;^.)j^M)ii)p,...PMi3) 1- 

x,,...p,,(l| X.... (3|5f)j/iM+i)...5f.^J/iL)|l)p,...pM|2)^-- (139) 

/Pl-'-PA/ TP1---PM 

Xpi...pA,(l| 1...L {2\Ti,„L\l)q^...qj,i\2) 1...L . 

^pi---pm ^n---iM 

We now decrease the number of relevant vector spaces in each of the above pseudo basic DWPF 
expressions. 

6.1.1 Reducing the number of relevant vector spaces - I 

Beginning with p^...p^^(l I {3\b[^\{ij,m+i) ■ ■ ■ b[^\{ij-l)\'^)p^...p^i\2) , we use elemen- 

tary matrix algebra and the fact that p(l|ep = to discard the values pi, . . . ,pm in the sum- 
mations of the ^{fi) operators to obtain, 

~ (2) ~ (2) 

Pi...PAf(l| ^ ^-^ (3|5i...L(AiAf+i) • • •^i...l(/^l)|1)pi...pa./|2)^ 1-i 
= nf=M+injii^3i(/ij,^pJpi...PA./(l| {3\B^^,„L (/iAf+i) • • • 5 (^l)|1)p,...p„|2) 

^P1---PM #Pl-..PAf ^Pl---PM 7^P1---PM 



(140) 



' (2) 

Since the operators B i _c (fi) act trivially on the vector space l^., i = 1, . . . , M we decrease the 
number of relevant vector spaces in the reference states by M - hence Eg. (11401) becomes, 

Pl-PAf(l| (3|^f x(^A/+l) • . .S|^.^.i(/iL)|l)pi...pAl|2) l...i 

7^Pl•••PA/ "FPl-'-PM 

= ntM+in-li&3i(/i.,ep,) ifiM+i) . . . B^\.., (^l)|2) 1.... 

^Pl---PM T^Pl-.-PA/ J^Pl---PM *Pl---PM 



I{i=M+l ni=1^3l(/^i,^p,)2'^_'M({/^fc}fc=A/+l,...,L|{6} J=i. -.-C' )• 
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6.1.2 Reducing the number of relevant vector spaces - II 

We now focus on the term i,„L{5\B^i\{fii) ■ ■ ■ b[^\{ij.m)\^)pi...pm\'^) ^ • We reiterate that 
due to the elementary relations in Eq.f llOip . there is no possibility of the Weyl matrices e^^ (^e^^ 
in the expressions for the operators producing anything but zero, hence they can be discarded 
from the calculations. Additionally, using the fact that ep'^^^|3)p = 0, we can discard the values 
I ^ Pi, . . . ,Pm iTL the summations of the B^lj^ operators. Applying the above facts, and elementary 
matrix multiplication, the expression now becomes, 

i...L(3|-B|^.^i(/ii) . . .s|^^^(/iM)|l)pi...pM|3)^ i- i 

X i...L(3|5pi!^p\,^(/ii) . . . Bp^^.^plj (/iM)|l)pi...pA/|3) 1...L , 

where Bp\'^^Jj^j{fi) is given by Eq. (11261) . 

Since the product nf=i 11^ , — ft is independent of the value of h, . . . ,1m, 

7^Pl,...,pjYj 3ll^j!^p^^j 3VSP;^5SjJ 

and the operators 5pif.p\^(/i) act trivially on the vector space Vi, i ^ pi, . . . ,Pm, we can decrease 
the number of relevant vector spaces in the reference states by L — M and perform the following 
simplification to Eq. fll4ip . 

i...^(3|5«J/xi) . . . 5«J^m)|1)p....p„|3) 1.... 

9^P1---PM 

~ llj=lll^ ^ = 1 b-i-,(£i £„ )e-,(£„ £,) Z^'1'--''m = 1 pi...pAl\<J|-Dpi---PA/lPlJ ■ • • -Dpi...pM lPMj|-'-/pi...pA/ 



[£l ^^^p,'ZlpM ''3lfeXtmL,5.) Pl-PAf(3|5pi^PA/(/^l) • • • ^m^..PAf(^Af)|l)pi...pM 



Applying all the results from this section we obtain the following form for the mixed DWPF of 
type B completely in terms of Boltzmann weights and the J-'-matrix sandwiched between reference 
states, 

^i^M({/^}> {O) = Zll<pi<...<pM<L Pi...PAf(lL 1 (2|J^l...L|l)gi...gAf|2)^i-i 



'(■B-l)a„ 1 ^ ..lit 1 . ..W(^'2), 



^^m' ({Aim}m=l,...,M|{'Cp„}n=l,...,A/)^LJAf({/^r}r=M+l,..,L|{'Cs} J'=l>- ->i' )■ 

#Pl.---.PAf 

(142) 
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6.2 Mixed DWPF for C 

We now apply the J-'-matrix to Eq. fll32p to twist the monodromy operators to obtain, 

(143) 

Similarly to the previous mixed DWPF we note that there is an J^i^^ operator on the far left 
of Eq.f ll43p as there is no equivalent expression such as Eq.f ll02p for mixed transpose reference 
states. As in the previous calculation for Z^j^lj we proceed by inserting two complete sets of states 
into Eq.f ll43p in separate stages for the sake of clarity. Consider first placing a complete set of 
states in between the operators C'l.xI^m+i) and C[^Il{i^m) to obtain, 

Z^lUMAO) = Ei<.,<...<..<lpi...pm(1| (3|C'?.^JM---C'^.^l(^i)|3)i...l 



/pi •••PA/ 

x,,...,,,(l| .... (2|J--i^Cf)^(z/z.)...Cf.)^(z/M+i)|l)p,...p,,|3) ..... . 

(144) 

To verify that the above terms of the complete set of states are the only ones which are non zero, 
consider the expression, 

C';'.^l(^m)...C'«,(^i)|3)i...l, 

and recall from Section [5731 that only the Weyl matrix (as opposed to e^^ (g) ej^) gives non zero 
terms when applied to the reference state |3)i...l- Hence since eP|3)i = |1);, this means the only 
terms in the complete set of states which are non zero are, 

z2 |1)pi-.pm|3)_^ i-i pi...pm{M ^^-^ (3|- 

l<Pi<-<Pi\i<L 

We now insert another complete set of states in between J-'i^^ and C[^\{i'l) in Eq. fll44p to obtain. 



^ P\...PM ' 



1| ..... (3|C|^/Jz/M)...q:/J^i)|3)i...L (145) 
x,,...,,,(l| ..... (2|C'f.)jz/^)...(7;;^)jz/M+i)|l)p....p„|3) ..... , 

i^Pl-'-PM "FPl-'-PM 

and proceed to decrease the number of relevant vector spaces in each of the above pseudo basic 
DWPF expressions. 
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6.2.1 Reducing the number of relevant vector spaces - I 

Beginning with I {2\Cf\{vL) . . .Cf\{vM+i)\l)p^,,,p^,^\'i) i...l , we use the fact 

^Pl---PM i=Pl---VM 

that Cp = to discard the values (pi, . . . ,Pm) in the summations of the C\ operators. 
Hence applying this fact, and elementary matrix multiplication, the expression simplifies to, 

= riilM+l rijll ^21 Pi. .. PM (1 L (2|C'*-^l..L [vl) ■ ■ ■C''\..L ('^M+i)|1)pi...pm|3)^ • 

(146) 

Since the operators C i...^ (z/) act trivially on the vector space Vp^, i = 1, . . . , M, we can decrease 

^Pl---PM 

the number of relevant vector spaces in the reference states by M - hence Eg. (11461) becomes, 

Pi-Pm{M (2|C'f.x(^i) • • •C'L.x('^Af+l)|l)pi...pM|3)^ ^-^ 

^Pl---PM ^P1---PM 

= nf=AmnSi&2i(i^.,ep,) (21^^.... k)...^^.... (z/m+ois) .... 

^P-L---PM ^Pl---PM 7^Pl---PM *P\---PM 

= ]\^=M+l ntil b2l{Vi,ip,)Zfl%{{Vk}k=M+l,...A{il} )■ 

6.2.2 Reducing the number of relevant vector spaces - II 

We now focus on the term p^...pj^ (1| 1...L {3\C[^\{um) . ..C[^\{iyi)\3)i,„L. Recall that due to 
the elementary relations in Eq. dlOip . there is no possibility of the Weyl matrices ej^ (g) ej^ in 
the expressions for the C[^\ operators to produce anything but zero, hence they can be discarded 
from the calculations. Additionally, using the fact that p{3\e'^^^ = 0, we can discard the values 
I ^ Pi, . . . ,pm in the summations of the C[^l]^ operators. Applying the above facts, and elementary 



matrix multiplication, the expression now becomes. 



p,...p,,(l| ..... {3\C[']^iu,j)...C['!M\3)^...L 

^Pl---PM 

Et...M^i n£i a3(^.,o)^3(o>ep,) (147) 



^Pl-PAliM I'-C- (SIC^Pi'^Pm^Z^m) • • • <^pl!^PM('^l)|3)l...i 



7^P1---PM 



where Cpi%{j {i')is given by Eq. fll20p . 
and the operators C'pi^pm('^) ^-ct trivially on the vector spaces Vi, i ^ pi, . . . ,Pm, we can decrease 



Since the product YlfLi Yl^ •'=1 (^3(^1^ ^j)^3(^j) ^pi-) is independent of the value of /i, . . . , Im, 
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the number of relevant vector spaces in the reference states hj L — M and simphfy Eq. fll47p as 
follows, 

Pl-PM 

= ni!!l •'^^ '^3('^i,'Ci)^3('Ci;^pJpi...pM(l|C'pi^..PM('^Af) • • •(5pi^Pm('^i)|3)pi...^^^ 

#Pi,...,PAf 

= n,= int a3('^i,^i)^3(^i,^pJ^M'^H{'^i}i=l,...,A/|feJj=l,...,Af)- 

Hence we obtain the following form for the mixed DWPF of type C completely in terms of 
Boltzmann weights and the inverse J-'-matrix sandwiched between reference states, 

Z^m{{'^},{^}) = Z]l<pi<...<p„<L 9i...<?A/(l| 1 (2|^L.x|l)pi...PAf|2) i- i 



T^Pl PM 

X^m'^ {{^m}m=l,...,M\{^Pn}n=l,...,M)^L-Mi{^r}r=M+l,...,L\{^s} s=i,...,L 



(148) 



7 Conclusions 

In this article we have argued that the factorized F-matrices associated to the i?-matrix of 
the U{1)^'^~^^ vertex model can be constructed for arbitrary Boltzmann weights. Our analysis 
is purely algebraic relying only on the structure of the i?-matrix as well as on the corresponding 
unitarity and Yang-Baxter relations. We have applied this formulation to the = 3 fifteen-vertex 
model which is the simplest extension of the asymmetric N = 2 six- vertex model in such family of 
integrable models. For = 3 we have exhibited the algebraic expressions of relevant monodromy 
matrix elements in the F-basis. This allowed us to compute the domain wall partition functions 
related to the creation and annihilation fields for arbitrary weights. 

We have motived our approach by showing that the Yang-Baxter relations for = 3 hide a 
general structure of Boltzmann weights. The underlying algebraic variety is at least governed by 
the intersection of two quadrics in the projective P'^ space leading us to a surface of del Pezzo type. 
Interestingly enough, this type of variety also governs the integrability of the N = 2 vertex model. 
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In fact, one can show that the weights of the asymmetric six- vertex model he on a cubic del Pezzo 
surface. A natural question to investigate is whether the del Pezzo structure persists for general 
N or even higher dimension algebraic variety emerges when > 3. In any case, this observation 
emphasizes the importance of attempts to establish results for integrable vertex models that are 
independent of any specific parameterization of Boltzmann weights. 
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Appendix A: The Yang-Baxter for = 3 

Here we describe some details entering the general solution of the Yang-Baxter equation for 
A^ = 3 exhibited in Section O The algebraic solution consists in the elimination of the weights 
dependent on the variables .^i and ^2 leading to determine the algebraic invariants constraining the 
remaining Boltzmann weights. We start by solving the relations involving only two triple prod- 
ucts, see Eqs.dmini). After eliminating the weights C2i(6,6), £32(^1, 6), &3i(6,6), ^13(6, 6), 
013(^1,^2) and &23(^i,G) one finds that there exists only three independent relations. They are 
easily separable providing us the following invariants, 

TTrT\=°^' h (c c ^ = '^2, , . N = '^3 tor 2 = 1,2, A.l 

&12(6,6) &2l(6,6) &23(^i,6) 

where ^i, 82 and ^3 are free parameters. 

Taking into account this solution the number of relations with three triple products ( IT41I2T1) 
reduce to thirty independent functional equations. Among of them there exists eight relations 
which are suitable for carrying out further elimination of weights. Their explicit forms are, 

012(^1,6)01(^1, ^3)^21 (6,^3) + &2l(6,6)Cl2(6,6)c2l(6,6) 

= ci2(ei, 6)^21(^1, e3)«i(6, 6), (A.2) 
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fc2l(6,6)Cl2(6, 6)^12(6, 6) + Ci2(6,6)ai(^l,6)Cl2(6,6) 

= ai(^i,6)ci2(6,e3)ai(6,6), (A.3) 



&12(6,6)ai(6,6)c2l(6,6) + C12(^1,6)C21(6,^3)&12(6,^3) 

= «l(ei,6)&12(ei, 6)^21(6, 6), (A.4) 

012(^1,6)02(6,6)^12(6,6) + &i2(6,6)ci2(6,6)c2i(6,6) 

= ci2(6,6)&i2(6,6)«2(6,6), (A.5) 

C23(6, 6)«3(6, 6)^21(6, 6)^23(6, 6) + &2l(6, 6)C23(6, 6)^23(6, 6)C32(6, 6) 

= 023(6,6)^23(6, 6)^21(6, 6)«3(6, 6), (A.6) 

'^l&2l(6, 6)C23(6, 6)^23(6, 6)^12(6, 6) + C23(6, 6)«3(6, 6)&2l(6, 6)C23(6, 6) 

= 03(6, 6)023(6,6)^21(6,6)^3(6, 6), (A.7) 

(^1023(6, 6)C32(6, 6)^12(6, 6) + ^32(6, 6)^3(6, 6)C32(6, 6) 

= 5ia3(6, 6)^12(6,6)032(6, 6), (A.8) 

(^1023(6, 6)^2(6, 6)^12(6, 6) + ^32(6, 6)C23(6, 6)C32(6, 6) 

= 5iC23(6,6)&i2(6,6)a2(6,6)- (A.9) 

From Eqs. (lA.2IIA.4l) we are able to eliminate the weights &2i(6,6), (^1(6,6), ^12(6, 6) and 
by substituting them in Eq.f lA.SP we find that it leads to the following constraint, 

ai(6,6)«2(6,6) + ^12(6,6)^21(6,6) -ci2(6,6)c2i(6,6) . . • 10 1n^ 

= ^4 for 2 = 1,2, (A.IO) 

ai(6, 6)012(6, 6) 

where 64 is a constant. 

The same procedure can be implemented for Eqs. (]A.6tlA.9l) . By eliminating the weights 
03(6,6), ^32(6,6) and 023(6,6) one obtains the additional constraint, 

«3(6,6)«2(6,6) + '5i&i2(6, 6)^23(6, 6) - C23(6,6)c32(6,6) 



03(6,6)^12(6,6) 
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55 for i = 1,2, (A.ll) 



where S5 is a free parameter. 

At this point we are left to eliminate only the Boltzmann weights 013(^1,^2) and 02(^1,^2)- 
This is done by performing linear combinations among certain remaining three terms relations 
coming from Eqs. (11411211) . Remarkably enough all the consistency conditions are solved by means 
of the following extra invariants, 

, c\h tc c\ = "^6 for 1 = 1,2, (A.12) 

&23(6,6)&2i(6,6) 

'^i [ai(6, 6)c23(6, 6) - ci3(6, 6)c2i(6, 6)] ^ 

C23(6,6)&2l(6,6) 

5iQi(6, 6)ci2(6, 6)^23(6, 6) - ci3(6, 6)^21(6, 6)c32(6, 6) ^ 
ci2(6, 6)^23(6, 6)^21(6, 6) 

5iai(6, 6)ci2(6, 6)^23(6, 6) - ci3(6, 6)&2i(6, 6)c32(6, 6) 



J7, 



(A.13) 



^8, (A.14) 



(A.15) 



ci2(6, 6)^23(6, 6)^21 (6, 6) 
Si [qi(6,6)c23(6,6) -ci3(6,6)c2i(6,6)] 
023(6,6)^21(6,6) 

^1 [qi(6, 6)032(6,6) - ci2(6,6)c3i(6,6)] ^ 

(^21 (6,6)c32(6,6) 
5iai(6, 6)^23(6, 6)c2i(6, 6) - C23(6, 6)^21(6, 6)c3i(6, 6) 
^23(6, 6)^21 (6, 6)c2i(6, 6) 

5iQi(6, 6)^23(6, 6)c2i(6, 6) - C23(6, 6)^21(6, 6)c3i(6, 6) ^ 

&23(6, 6)&2l(6, 6)C21(6, 6) 

Si [ai(6,6)c32(6,6) -ci2(6,6)c3i(6,6)] _ . 

621(6,6)032(6,6) - ^^-''^ 

where ^e, (Jy, ^g, Sg and 5io are yet new free parameters. 

It turns out that the remaining Yang-Baxter functional relations lead us to branches that 
impose further constraints among certain weights and the invariants obtained so far. We find that 
one such possible branch is, 

^2(6, 6)621(6, 6) _r Q2(6, 6)623(6, 6) _ 61(6467-61) -^lO 17^ 

ai(6, 6)612(6, 6) '«3(6, 6)612(6, 6) 5? ' ' ' ^ ' 
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while the invariants S^^S^^Sq^Ss^Sq^ Siq are fixed by, 

^3 = r /r / — FT' ^5 = '^4, 4 = 0, Ss = St, 6g = 6io = — —. (A.18) 

02(^0407 — 01) O4O7 — Oi 

We have now reached a point where all the weights in the variables ^1 and ^2 have been 
eliminated, while the weights in the variables ^1,2 and ^3 are constrained by the algebraic invariants 
( 1A.lllA.10tlA.17l) . The final step of our analysis consists to make the intersection of these algebraic 
invariants. The procedure for performing such intersection is as follows. We first note that the 
weights a2(^i, 6), &23(^i, ^3), &3i(6, 6), ^32(^1, 6), 03(6, ^3), £13(6, 6) and 032(6, 6) can be linearly 
extracted from flA.l|A.imA.17p . leading us to, 

^ijSi - '^4'^7)Qi(6,6)ci2te,6)c2ite,6) (AiQ) 

""'^^''^'^ [5,aite, e3)-M2ite, 6)] [(5i-Waite,e3) + 57^21^,^3)]' ^ ' ^ 

, (f f \ ^ (^1 - ^4^7)^21(6, 6)c23(6,6)c3lfe, 6) 

''^^"^'^ C2ite,e3)5l[(5l-5457)aite,e3) + '57&2lte,e3)]' ^ ' ^ 

&3l(6,e3) = 52&2lte,e3), (A.21) 
, (<: )^ '^l ^7^21 (6, 6)Cl2te, 6)c2lfe, 6) . . 291 

^^'^^^'^^^^ [5,aite,e3)-M2ite,e3)][(-5i + 5457)aite,e3)-M2i(6,e3)]' ^ • ^ 

^ X - (54^7)Qlte,6)c23te,6)c3l(6,6) , . 2-?^ 

""'^^"^'^ C2ite,e3)[(<5l- Wai(6,e3)+M2lte,e3)]' ^ • ^ 

C'lC2l(4i,^3) 

^ ^ _ ('^1 - ^4^7)012(6, '^3)c3ite, 6) /A 

where i = 1,2. 

We next substitute the weight (1A.19I) in Eg. ( lA.lOl) and as result we find that the weights 
aite,^), &i2te,6), &2ite,6), ci2(6,6) and C2i(6,6) are constrained by, 
Si{5i57 - 5i] 



5? 



^12(6, 6) - (^401(6, 6)^12(6, 6)&2ite, 6 



+&12(6,e3)&2lte,e3) -&2l(6,6)Cl2(6,e3)c2lte,e3)=0, for Z = 1, 2. (A.26) 



By performing the definition ^7 = 61 Ai and ^4 = A2 we see that the form of Eg. flA.26p is the 
same as that of the hypersurface (1251) given in the main text. This is the case because the spectral 

47 



parameter ^3 is a common variable for all the weights entering Eg. (1A.26|) . Therefore, through the 
identification, 



02^,6) = ate), &12te,6) = K^i), &2l(6,6) = K^i), Ci2te,6) = Cte), C2lte,6) = Cte), 
C23te,e3) = /^lte), 031(6,^3) = /^2te), (A.27) 

we see that Eq. flA.26P becomes exactly Eq. fl23l) . 

We conclude by observing that the earlier U{1) U{1) Yang-Baxter solution presented in the 
literature [201121] is indeed a particular case of the i?— matrix given in the text. In fact, the so 
called Perk-Schultz solution associated to the Ug[Su{3)] quantum algebra is obtained by setting, 

5^=82 = 1, Ai = g, A, = q + l/q, a(0 = l, = = 9^^, 

r - r 

hiO = ^h^iO = c(0 = c(0 = e^^, (A.28) 

where g is a free constant and ^ is the spectral parameter. In this special case we see that the 
i?— matrix is of the difference form, i.e., R = R{C,i/^2)- 



Appendix B: The analysis for L = 2 

Here we present the explicit expressions of the entries fl91f92l) together with the corresponding 
simplifications using the Yang-Baxter and unitarity relations; 

C2lte,6) {^3l(At,^l)&32(Ai,6) - &32(yU,^l)&3l(yU,6)} = 

V ' 

apply Eg. 11211-13.1.21 
&3l(At, 6)a3(At, 6)c3lte, 6) + Cl3(/i, 6)c3l(/U, 6)&3lte, 6) 

" V ' 

apply EQ.lfT4ll-l3.1l 

-a3(/U,6)&3i(/i,6)c3i(6,6) = 

l^f^ = ^32(/i, 6)«3(At, 6)c32te, 6) + C23(At, 6)c32(/i, 6)^32(^1, 6) 

^ V ' 

apply Eq.|[T4ll-l3.2l 
-a3(/i,^l)632(/i,6)c32(6,6) = 
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&2l(Ai,6)c32(Ai,6) 

apply Eq.tT0l--f3.2> 

&32($2,gi) fc32(/^, 6)Q2(/^, ^i)&32(6, 6) + &32(/^, ■Cl)c32(/i, 'Cl)c23(6, Cl)} 

apply Eq.lfT9t-l3.2j 

C32 (/^,g2 )fe32 )a2 ($2 ,gl ) 

C32(Ai, 6)a3(yU, 6)C32(^1, 6) + &23(yU, 6)c32(/i, 6)^32(6, 6) 

^ . ' 

apply Ea.lfT6ll-l3.2l 

«3(/U, 6)c32(/i, 6)03(6, 6) 

C32 (^.^1)^31 )fc2i (gi ,6) 

b3i{a,6) 

C32(/^,gl)fc32(/^,6) 
^32(51,52) 

C32(/i,6)a3(/i,6) 

&2l(/i, 6)c32(/i, 6)C21(6, 6) + Cl2(/^, 6)c3l(/i, 6)^21 (6, 6) 

apply Eq.JT3j-{2,3,l} 
, ^ s 

c32(a',6)&3i(/^, 6)^21(6, 6) C3i(6.6) 
631 (?i ,6) 

&2l(/i, 6)C32(/U, 6)C21(6, 6) + Cl2(/i, 6)c3l(/i, 6)^21 (6, 6) 

^ V ' 

apply Eq.|[2T1l-i3.2.1j 
-C32(/i,6)&2l(/i,6)c3l(6,6) = 
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&3l(/^,6)c23(/^,6) 
&32(/i,6)c23(/^,6) 

«3(M.gl)c23(M.g2) 
b32(€l,$2) 

C23 )fe21 (At,g2)fa31 (€l ,6) 

^21 (6 ,6) 

C23(/i, 6)«2(/i, 6)^32(6, 6) + &32(/i, 6)c23(/i, 6)C32(6, 6) 

V ' 

apply EQ.lfT8ll-l2.3l 

C23(/^, 6)&32(/i, 6)«2(6, 6) 

apply Eq.tTni--f2.3> 

C23(/^,ei)&23(/i,6) - °3(M.6)c23(M.6)^^C32(ei,6) 

b32(g2,6) i^23(/i, 6)c23(/i, 6)^32(6, ^1) + C23(/i, 6)Q3(/i, ^l)c23(6, 6)} 

apply Eq.|[T7t-l2.3) 

03 )C23 )a3 ($2 ,gi ) 
fe32(6>5i) 

Cl3(/i, 6)C21 (yU, 6)^31 (6 , 6) + &31 (yU, 6)c23(yU, 6)C31 (6 , 6) 

apply Eq.tT3l-(3.2.1> 



C23(M,€l) hlif^, 6)^31(6. 6) C2l(gl.6) 

&2i(6,6) 

Cl3(/i, 6)c2l(/U, 6)&3l(6, 6) + hlifJ', 6)c23(/i, 6)c3l(6, 6) 

^ V 

apply Eq.|[2TJ-{2,3,l} 
-C23(Ai,{l)&3l(Ai,6)c2l(6,6) = 
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apply Eq.tT0l-(3.1> 

apply Eq.lfT9ll-l3.1l 

C31 (^.6)^31 (fi,£.i )ai ($2 
fe3l(6,?l) 

C2I (/i, ^1)C32 (/i, 6)C21 (6 , 6) + ^12 (/i, ^1)C31 (/i, 6)^21 (^1, 6) 

apply Eq.tT(n--f3.2> 

C3l(/^,6)fe32(/^,g2)c2l(gl.6) £32(^1, ^2) 
632(6.6) 

g^t|^{c32(/X,e2)c2l(/i,ei)fe32(6,ei)+fe32(/i,6)c3l(/^,^l)c23(6,^l)} 
^ V 

apply Ea.|[20ll-i3.2.1j 

+&12(/i, 6)c3l(/i, 6)&2l(6, 6) 

g|f^{c2i(6^6)c^^^ 

apply Eq.^-{2,1} 

g||^{&32(/^, + &2l(6, 6) [&12(/^,ei)&32(6,ei) - &32(/i,ei)&12(6,ei)] } 

apply Ea.lfnt-il.3.2j 

fe32(/^,6)c3l(/^,6) 
^32 (6 .6) 

C3I (/i, ^l)a3(/U, 6)C31 (6 , 6) + &13 (/i, ^1)C31 (/i, 6)^31 (^1, 6) 

" V ' 

apply Eq.|[T6ll-i3.1> 

a3(/^,^i)c3i(/i, 6)03(6, 6) 

C3l(M,6)fe3l(/J,6) 

^31 (6 ,6) 

£31(^.6 )fe32(/J,6) 
f'32(6>6) 

C3l(/i,6)«3(/^,6) 

apply Eq.(Tn}-{3,2} 

C2l(/i,6)c32(/.,6) - -3.(^-^-)^32(.6)^^^C32(6,6) 

fe32(6,6) 6)c2l(/i, 6)^32(6, 6) + ^32(/i, 6)c3l(/i, 6)^23(6, 6)} 

apply Eg. 1201-13. 2.1> 

C31 (/^,6 )fe32 )C21 (6 .gi) 

632(6.6) 
C31 (M.^I )fe32 (/J,6 )C21 (gl ,6 ) 

632(6.6) 
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Z^!^'^ = &3l(/^,^l)Cl3(/i,6) 

^(-Bl) _ b32(M,€l)ci3(M,g2) 

- b2l(6,6) 

^(-Bi) _ a3(/J,gl)ci3(A'.6) 
'^S b3i(Ci,6) 



Cl3(/i, 6)^31(6, 6) + &3l(/U, ^l)ci3(/i, 6)c3l(6, 6 



apply Eq. |[T8)| -I1.3> 



= Ci3(/i,^i)63i(/X, ^2)01(6, 6) 

4^'^ = - T^fc'i^l {c23(/^, 6)Cl2(/^, 6)^32(^1, 6) +&32(/i,^l)Cl3(/^,6)c32(6, 6)} 

V ■' apply EQ.|[2T1l-il.3.2j 

apply Eq.|[T0t-l2.1) 

+Cl3(/i, 6)&12(/i, 6)^32(6, 6) 

= gfe|\ { Cl2(ei, 6)C12(6, 6)^ &32(;U, 6) + fel2(/i, 6)^32(6, 6)^21(6, ^l)} 
apply Eq.ll9ll-{1,2} 

= afelr i^32(/^, 6) + ^21 (6, 6) pl2(y^, 6)^32(6, 6) - &32(/X, 6)&12(6, 6)] } 

apply Eq.|[T3t-n.3.2> 

_ Cl2(/J,6)fc32(At,g2) 
b2l(6,6) 

apply Eq.tT0l--f3.1> 

JBl) _ „ (,, t \h (,, t \ °3(M.gi)ci3(/^.g2) £31(^1,^2) 



= baife.gi) {bisif^, 6)ci3(/i, 6)^31(6, -^l) + Ci3(/i, 6)Q3(/^, 6)^13(6, 6)} 

apply Eq.llT7t-il.3) 

_ ^3 (/^,g2 )ci3 (m.^i )a3 (6 ,€1 ) 

631 (6, a) 

apply Eq.tT0l--f2.1> 

^(Bl) _ fe32(M,gl)ci3(M2) C2l(Cl,6) 

7 b2i(a,52) 

_ &23 (A'.a )C13 )C12 (6 ,gl ) 

^21 (6,?!) 

4^'^ = b2i(gi.g2) {^23(/i, '^l)ci2(/i, 6)^32(6, 6) + ^32(/i, ■Cl)ci3(/^, 6)^32(6, 6)} 

apply Eq.|(2lj-{1,3,2} 

_ Ci3 )fe32 (At,6)ci2 (6 ,6) 

^21 (6 ,6) 
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Appendix C: The analysis for L = 3 



In what follows we present the explicit expressions of the entries ( l93l) together with the corre- 
sponding simplifications using the Yang-Baxter and unitarity relations; 

apply Eq.lT3j-{3,2,l} apply Eg. tTol--f 3.21 

(Ci) ^ C2l(M,a)c32(/^,6) ^31 if^, ^3)^21 (6, 6) _ C3l(/^,a)fe32(M2)b3l(M,g3) £32(^1, 6) 

"1 f'2i(Ci,€3)fe3i(6,6) 631 (Ci ,6)632 (€1,6) 

apply Eq.lT3)-{3,2,l} 



l{^i,il)c32{^l,a)b2l(p,i3) I C31 (A'.gl )b32 ) ^31 (jU, 'Cs) ^21 (^1 , C23 (6 .gl ) 

621(6,6) 632(6,6)631(6,6)621(6,6) 



621 (m,6) 



{c32(/i, 6)c2l(/i, ^l)fe32(6, ^1) + &32(/U, 6)c3l(/i, ^l)c23(6, 6)} 



632(6,6)621(6,6) 

apply Eq.|[2011-I1.3.2) 

_ 621 (m,6)c31 (m,6 )632 (At,6 )C21 (6 ,■^1) 

632(6,6)621(6,6) 

apply Ea.tT0l--f3.2l 

^{Cl) ^ C2i(/J,6)c32(m,6)632(/^,6) _|_ C3i(M,6)632(m,6)632(m,6) £32(^1, ^2) 
2 632(6,6) 632(6,6)632(6,6) 



632(^,6) 



{c32(/i, 6)c2l(/i, 6)fe32(6, ^1) + &32(/U, 6)c3l(/i, 6)C23(6, 6)} 



632(6,6)632(6,6) vi 

apply Ea.|[20ll-ll.3.2) 

_ 632 (m,6)c31 (A',6 )632 (At,6 )C21 (6 ,€1) 

632(6,6)632(6,6) 

apply Eq.tT0l-j"3.2> 

= C2l(/.,6)c32(/i,6)a3(/.,e3) - -^-(^'^^)^32(M2)a3(^^^^^^^^^^ 032(^1,6) 

= ^t||$y {c32(/^,6)c2l(/i,6)fe32(6,6)+fe32(/i,6)c3l(/i,G)c23(6,ei)} 

apply Eq.CT-n.3.2> 

_ Q3(A',6)c3i(m,6)632(m,6)c2i(6,6) 
632(6,6) 

apply Eq.(l3)-{3,2,1} apply Eq.tTol-f 3.2> 



o{Ci) ^ C2i(m,6)62i(m,6)c32(m,6) _ C3i(M,6)fe3l(At, 6)^2l(6^ 6) 632(^^,6) C32(G, 6) 

^1 621(6,6) 631(6,6)621(6,6)632(6,6) 

= &2i(6'6fe?(6,6) fc32(/i, ^3)c2l(/i, 6)fe32fe, 6) + fe32(/^, 6)c3l(/i, 6)c23fe, 6)} 

apply Eq.ll20t-n.3.2> 

_ 621 (/j,6 )c3i (^,6)632 (/j,6 )c2i (6 ,6) 
621(6,6)632(6,6) 
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apply Eq.tTot-f 3.2> 



apply Eq.JT0}-{3,2} 

C3l(M.gl)&32(M,g2)fe32(M,g3)c32(6.g2) £32(^2, ^3) C31 )fe32 (/J,6)fe32 (At.g3)a2 (gl .6)c32 (gl .6) 



^32 (5l ,6)''32 (6 ,6 ) f'32 (6 ,6 )ft32 (?1 ,6 ) 

ggt) i^32(/i, 6)a2(/i, 6)&32fe, 6) + h2{fi, e3)c32(/i, 6)c23(e3, 6)} 

apply Eg. 1191-13.21 

~ blt\ciil)b32{S3^^^^^ ) i ('^1 ' ^3 ) C32 ('Cl ; 6 ) C23 (6 ; 6 ) + C32 (6 ; 6 ) Q2 (6 ; 6 ) ^32 ('^3 , 6 ) } 

apply Eq.llT9t-l3.2j 

apply Ea.tT0l--f3.2> 

C21 (M.gl)c32 (^.6)^32 (M.6)a2 (6,6) _ C3l(/^.6)fe32(M.6)fe32(/J,6) C32(^l,^3) a2(6,6) 
b32(6>6) b32(5l,C3)t'32(6i6) 

= 63:fe:|)£i"g) i"32(^, 6)021 (/i, 6)^32(6, + fe32(/i, e3)c3l(/X, ^1)023(6, ^l)} 

apply Eg. 12011-11.3.2) 

fe32(/J,6)a2(6.6)c3l(M.6)fe32(/J,gl)c2l(6.6) 
b32(6i6)t'32(6i?l) 

C2l(/i, 6) {C32(/^, 6)a3(/i, 6)C32(6, 6) + &23(/i, 6)c32(/i, 6)^32(6, 6)} 

apply EQ.lT6t-i3.2) 

- '"^"^^Crit'^'^ i^32(^, 6)a3(^, 6)032(6, 6) + C23(/i, 6)c32(/i, e3)&32(6, 6)} 

apply Eg.lTH-i3.2) 

apply EQ.tT0l-(3.21 

I C3l(M,a)a3(M.6)fc32(/^,6) £32(6,6) C32(6.6) _ C31 )a3 )fe32 (m.6)c32 (gl ,6)^3 (6 .6) 

632(6,6) 632(6,6) 



C2l(/^, 6)a3(/^, 6)c32(/i, 6)^3(6, 6) - 



C31 (M,gl )C32 (gl ,6)«3 (^,6)632 (m,6)c32 (gl ,6) 



632(6,6) 

''-'^23(6, 6)032(6, 6)^32(6, 6) + ^32(6, 6)% (6, 6)032(6, 6)} 



V 

apply Eg.lTIt-l3.2> 

apply Eq.tT0l-l3.2> 

C2l(/i,6)a3(/i,6)c32(/., 6)^3(6, 6) - '''^'''^ C32(6,6)^^^^3(..6)632(.,6)C32(6,6) 

apply Eq.tT0t-(3.21 

C3l(/^,6)Q3(/J,6)632(M,6)a3(6,6)632(6,6) £32(6, 6) 
632(6,6)632(6,6) 

C2l(/i, 6)«3(/U, 6)c32(/i, 6)03(6, 6) 

i(/J,6)"3(m,6)632(/J: 
632(6,6)632(6,6) 



^"^t?.5!rd^l1.f {^32(6, 6)C23(6, 6)^32(6, 6) + &32(6, 6)a3(6, 6)C23(6, 6)} 

apply Ea.fT5ll-i3.2> 



54 



°'^fattf)'^'^ i^32(/X, ^s)c2li^, ei)&32fe, + &32(/^, 6)031 (/i, ^1)023(^8, ^l)} 

apply Ea.|[20ll-ll.3.2j 

"3(M2)"3(6,g3)c3l(/J,g3)fc32(/J,gl)c2ltei'?l) 
''32(?3,6) 

apply Eq.trnt-f 2.1> 



C2 l(M.gl)fe2l(/J,g2)c32(/J,g3) C2l(6, 6) 
b2l{6,€2) 



ai(Ai,6)c2i(/i,6)c32(/^,C3) - 

apply Eq.(in)-{3,1} apply Eq. (in)-{3,2} apply Eq . tTol--f 3.21- 

I C3i(At,6)b3i(M,6)fe32(M.€3) C3i(,^i,6) £32(6,6) _ ai(M,6)e3i(M,g2)fe32(At,g3) £32(6; 6) 
b3l(6,?2)fe32(a,6) 632 ,6) 

gf(|f\ {c2l(/X,6)ai(/i,6)fe2l(6,6)+fe2l(/^,6)c2l(/i,6)Cl2(6 

apply Eq.llT9t-i2.1) 

+£feifeiS{^3i(/i,6)c3i(/i,6)c3i(6,6) + c3i(/i,6)«i(/i,6)&3i(6,6)} 



apply Eq.lfT9ll-l3.1l 

apply Eq.tT3l-(3.2.l1- 



C32 (M.gs )C21 )fa21 )ai (6 ,gl ) I fc32(M,$3)c23(g3,g2)c3l(M,g2)ai(g2.gl) ^31 (/^; 6) 

b2l(6,5l) b32(e3,6)fe3l{6,6) 



{C32(yu, 6)c2l(Ai, 6)^32(6, 6) + &32(/i, 6)c3l(/i, 6)C23(6, 6)} 



fc2i(M.gi)ai(6.€i) 

f'2l(6:?l)''32(€3>6) 

apply Ea.|[20t-ll.3.2> 

fe2l(At,gl)ai(g2.6)c3l(A',g3)b32(At,g2)c2l(g3,g2) 
b2l(€2,Cl)fe32{C3,^2) 



&12(/i, 6)c2l(/i, 6)c32(/i, 6)C21(6, 6) 



apply Eq.tTol-f 3.2> 

fel2(At.gl)c3l(M,$2)fe32(M,g3)fe2l(gl,g2) ^32(6; 6) 
^32 (6 ,6) 

apply Eq.tT0l-l3.2> 



+C2l(/.,6)a2(/.,6)c32(/i,6)c2l(6,6) - "^^"'^^^"^^"'^^^^"^t?(Sf''^^ "'^^^^'^'^ 

apply Eq.tT0l--f3.2> 

I C3l(A'.€l)fe32(/^,g2)fc32(/J,g3)c32(6.g2)c2l(6,g2) £32(6; 6) _ C31 (At,$l )&32 ifify)b32 (A'.6)a2 (gl ,g2)c21 (gl .6)c32 (gl .gs) 
632 (6 ,6 )ft32 (6 ,6 ) ^32 (6 ,6 )fe32 (6 ,?3) 

'"^"^dSf/'^'^ i^32(/i, 6)C21(/X, 6)^32(6, 6) + h^if^, 6)C31(/X, 6)C23(6, 6)} 

apply Eq.|[20ll-il.3.2j 

+ '"^ fatVS^^^^^ {^32(/i, 6)«2(/i, 6)&32(6, 6) + &32(/i, 6)c32(/i, 6)c23(6, 6)} 

apply Eq.llT9t-i3.2) 

~ ^bstllu^iV^afi^^^^ ''^^ "* { ^32 (6 ; 6 ) C32 (6 ; 6 ) C23 (6 ; 6 ) + C32 (6 ; 6 ) Q-2 (6 ; 6 ) ^32 (6 ; 6 ) } 

apply Eq. 1191-13.21 
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bl2 (^,$1)621 (gl ,6)C31 (/^.g3)fe32 (m.6)c21 fo ,6) _j_ C21 (M.gl )C21 (gl ,6)^32 {^J.,^3)b32 (/J,6)a2 (gs ,6) 



621(6,6) b32(?3>6) 
apply Eq.tT0l--f3.2> 

C3l(M.gl)fe32(/J,6)fe32(/J,6)c2l(6.6) £32(^1, ^3) a2(g3.6) 
632(6:?3)632(C3:6) 
bl2 (^,6)^21 (gl ,6)C31 (M.6)fe32 )c21 (6 ,6) 

62l(6:C2) 
C21 (gl ,6)a2 ($3 ,6)^32 ) 



apply Ea.|[20ll-ll.3.2> 



apply Eq.tTSt-f 1.3.2> 



bujfJ', 6)^32(6, 6) fe2l(gl,g2)c3l(/^.g3)fe32(/^,6)c2l(g3,6) , C31 (/^.gg )fc32 )c21 (gl ,6 )n2 (6 ,6 )fe32 (^.gl )c21 (6 ,gl ) 

f'32(C3.gl)fe32{?3,6) 632 (6 .6)632 (6, 6) 



C31 (^.6)632 (^.6)632 (M.gl ) 
632(6:6)632(6:6) 



{^12(6, 6)C21(6, 6)&2l(6, 6) + C2l(6, 6)^2(6, 6)C21(6, 6)} 



apply EQ.lfT6ll-i2.1> 

£31 (^,6)632 (^.6)632 (^.6 )«2 (6 :6 )a2 (6 ,6)c2i (6 .6 ) 
621(6:6)632(6.6) 

apply Eq.tTol-f 3.2> 

_ 613 (M,gl)c21 (/J.6)C32 (m.6)631 (6 .6)632 (gl .6) 6i3(M,6)c3i(/^,6)632(/J:6)63i(6:6)632(6:6) 0^2(^2, 6) 



632(6:6) 632(6.6)632(6.6) 

apply Eq.tT0l-(3.2> 

I C31 (M.gl )C23 (/^.6 )C32 (^.6)^31 (gl .6 )632 (gl ,6 ) _ ^31 (M.6 )a3 (^.6)632 (/J.6)c31 (6 .6)632 (6 .6) ^32(6; 6) 
632(6.6) 632(6.6)632(6.6) 

I C31 (/J.gl )632 (A'.6)a3 (m.6)c21 (6 .6)C32 (gl .6) 
632(6.6) 

'"^'6;t\6^^^^^^^^^^^^ {^32(/X, e3)c2l(/i, 6)&32(e3, 6) + &32(/i, e3)c3l(/i, e2)c23(e3, 6)} 

apply Eq.|[20t-ll.3.2> 

+ '^"^"6;t\6^^^^^^^^^^^^ {^32(/i, e3)c23(/^, 6)^32(^3, 6) + &32(/i, e3)a3(/X, 6)023(^3, 6)} 

apply Ea.lfT5ll-i3.2> 

I C31 (At.6 )632 (m.6)"3 (^.6)^21 (6 .6)C32 (6 .6) 

632(6.6) 

613 (/^.6 )63i (gi ,6 )632 (6 .6)c3i (^.6)632 (m.6 )c2i (6 .6) 
632(6.6)632(6.6) 

"^'^'^^ 63f(6'!6)63f(6!i^^^^ {^32(67 6)^31 (6, G)C23(6, 6) + £32(6) ^3)^21(6; G)^32(6, 6)} 

apply Ea. l(20)l -ll.3.2> 

apply Eq.tT2t-(3.2.1t 

6i3(/^.6)632(m.6)c3i(M.6) ^3i(6;'C2) 632(6.6)c2l(6.6) _|_ C3i(/J,6)632(At.6)a3(M.6)c3l(6.6)c2l(6.6) 



632(6.6)632(6.6) 632(6.6) 

632(6.6) 



632(M.6)C21(6.6) {6,3(;x,ei)c3i(/i,e3)fe3l(ei,6) + C3l(/X,6)a3(/i, 6)031(^1, 6)} 



apply Eq.|[T51l-l3.11 

C32 (p.6 )c2i (6 .6 )«3 (^.6 )c3i (m.6 )q3 (gi ,6) 
632(6.6) 
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apply Eq.tT3l--(3.2.1> 



r(Cl) ^ C3l(Ml)b32(M,6)fe3l(/^, ^3)^21(6. ^3) C2l(gl,6) 
"1 b32{Cl,6)b3l{Cl,e3)fe2l(6,e3) 
_ C3l(At,gl)b32 (At,g2)fel (/J,g3)c21 (gl ,6) 
^32 (51,6)^21 (6 >&) 
^(Cl) _ C31 (At,gl )&32 (^.6)^32 (/^,g3)c21 (gl ,6) 

2 f'32(?l,?2)fe32(?l,6) 

^(Cl) _ C3l(M.gl)fe32(M.6)"3(M.g3)c2l(gl,g2) 

3 ^32 (6 .6) 

apply Eq.JT2}-{3,2,l} apply Eg . tT0l-<2.1 1 

j,{Ci) ^ C31 (M.gi )b3l (M,6)fe32 (gi .6)C21 (gi .6 ) _ C3i(M,gl)&32(/i,6)^3l(/i,6)c2i(6,6)fe2l(6.€3) 021(6,^3) 
^1 63l(?l,6)b32Kl,«3) f'3l(6,6)fe32(ei,6)t'2l(6,«3) 

= blt}cii2)bsl($u^3^^^ {^21(^1, ^3)Qi(^1) 6)^21(^3, 6) + ^21(^1, ^3)^21(6; 6)Cl2(^3; 6)} 

apply Eg. 1191-12. 11 

apply Ea.tT3t-f 3.2.1> 



C3I 



bsijfl, 6)^21(6. 6) fe32(/^.6)c2l(gl.g3)ai(g3.6) 



63l(6,6)fe32(5l,e3)b2l{C3,6) 
C31 )fe21 (At,g2 )b32 (/J,g3 )C21 (gl ,g3)ai (fe ,6 ) 

632(6:?3)''2l(C3,C2) 



(,{Cl) _ C3I (^,^1)^32 (^.6)^32 (M.g3)c21 (gl ,g3)a2 (6 ,^3) 



fe32(6:C2)ft32(5l,C3) 



^2 

_ C31 (M.gl )"3 (/J,g2)fe32 (/J,g3)c21 (gl ,6) 
3 fe32(6:C3) 

apply Eg.tTol-f 3.11- 

j{Cl) ^ ai (P,il )C31 )fa32 (/^,g3 )C21 (6 ,^3) _ C3l(M.gl)fe3l(M,g2)''32(/^.g3) C3l(,^l,^2) ^21(6.6) 

1 632(6:?3) 63l(?l,6)''32(6>?3) 

apply Eq.ll0j-{3,1} 

I C3l(At,gl)b32(M.6)fc3l(M.g3) C3l(^l,,^2) C2l(6.€3) _ C31 (/^,gl)fe32 (M.g2)''31 (M.g3)c31 (gl ,g3)c23 ($2 ,^3) 
631(6,6)632(6,6) 631(6,6)632(6,6) 
C2i(m,6)c32(m,6)63i(m,6)c2i(6,6)62i(6,6) 
621(6,6)631(6,6) 

= t^t^iS:!SS!) i"3l(/i, ^,)a,{f,, eO&3l(6, 6) + hlif^, 6)C31(/X, 6)C13(6, 6)} 

apply Ea.lfT9ll-l3.1l 

~ 631 (6 ,fo6t'2^(6^^^^^ {^13(^2, ^i)c2i(6) ^3)^31(^1; ^3) + ^31(6) ^1)^23(6) ^3)^31(6) ^3)} 

apply Eg. 12111-12.3. II 

C2l(M,6)c32(At,6)63l(A',6)c2l(6, 6)621 (6,6) 
621(6,6)631(6,6) 

_ 632 (At,6 )c2i (6 ,6 )c3i (m,6 )63i (/J,gi )ai (6 ,6 ) _ C31 (/j,6 )632 (^',6)631 (m,6)c32 (6 ,^^1)631 (6 ,6)^21 (gi ,6) 
632 (6 ,6 )63i (6 ,6 ) &31 (6 , 6 ) h2 (6 , 6) 631 (6 ,6 ) 

v ' 

apply Eq.tT2l--f3.1.2> 

C21 (M,gl )C32 (At,6)631 (^,6)^21 (6 ,6)621 (6 ,6) 
&2l(6,6)&3l(6,6) 

^ V ' 

apply Eq.tT3l-(2.3.1> 
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b32(A',g3)c2l(6.6)c3l(M.6)fe3l(/^.6)ai($2:6) 
632(^2, ?3)fe3l(C2,Cl) 

- ^:'tll£)^(^^!) {^32(^,6)c3l(/i,G)c23fe,ei) + C32(/X,6)C21(/^, 6)^32(6, 6)} 
^ V ^ 

apply Eq.|[20t-ll.3.2l 

apply Eq.tT^-f 3.1.2> 
. 



b32 (/^,6)c2l ($2 ,g3)c3i (^.6)631 (M.gi )ai (g2 ,gi ) _ C21 (6 .6) ^31 (A^; ^3)^32(/^i ^l) C3l(M2)c2lte.gl) 

632 (6 ,^3)631 (6 ) bsi (€1 ,6 )i'32 (6 ,€1 ) 

apply Eq.tT2l-(3.1.2> 



b32 (/^,6)C21 ($2 ,g3)c31 (Ai.^ )^31 (^.gi )ai (6 ,€l ) _ ^32 (^'.$3 )c3l )c2l (gl ,6)^31 (M.gl )c21 (g2 ,gl ) ^32 (^2 , '^3 ) ^31 (^2 , '^1 ) 

632 (6 ,^3 )fe31 (6 ,?1 ) b32 (6 ,6 )''31 (?1 ,6 )ft32 (6 ,51)^31 (6 ,5l ) 

apply Eq.JT3j-{2,3,l} 



fa32 (m,$3)c21 (6 ,g3)c31 )''31 (P'.gl)"! (6 ,gl ) _ ^32 (At,g3)c31 (M:6)fe31 )C21 (6 .g3)c21 (6 .gl) ^31 (^2 7 ^3)^21 (G, ^3) 

632 (6 ,6 )ft31 (6 ,Cl ) f'32 (6 ,6 )f'31 (?1 ,6 )^'32 (6 ,51)^1 (5l :53 ) 

apply Ea.tT0l--f2.1> 

fa32 (M,g3)c21 (6 ,g3)c31 (At,6)fc31 (/^.gl)"! (6 ,gl ) _ ^32 (At,g3)c31 (/^,g2)&31 (^,6 ) C2l{C,l, (,3) C21 ($2 .$1 )b21 (6 .6) 
632 (6 ,6 )^'31 (C2 ,Cl ) ^32 (C2 ,€3)631 (5l ,6 )f'21 (Cl ,^3 ) 

b2Ti^3il)b32(^2^^^ {^21 (6) ^3)0-1(6) 6)^21 (^3; 6) + ^21(67 ^3)^21(6) 6)Cl2(^3) 6)} 

apply Ea.llT9ll-l2.1j 

apply Eq.JT3j-{3,2,l} 
,. ^ ^ 

fe32(M3)c3l(M.6)&3l(At, ^1)^21(67 6) C2l(g2.g3)ai(g3.gl) 

b2lK3,6)fe32{6,6);'3l(6,«l) 
b32 (A',g3)c31 (M.g2)b21 (M.gl )C21 (6 ,6)"! (6 ,€l ) 
fel(C3,Cl)fe32(C2:C3) 

apply Eq.tT0l--f3.2> 

C21 (M.gl )C32 (At.g2 )fe32 (M.g3)c21 (gl ,g3)a2 (6 ,€3) _ C31 )fc32 (/J,6)fe32 (/J,€3) C32(^1,.C2) C21 (gl ,g3)a2 (6 .€3) 

f'32(?2:?3) 632(Cl,6)''32(S2,?3) 
I bi2 (At,gl)c3l (At,g2)fc32 (^.$3)^21 (gl ,g3)c21 (6 .gs) 
632(6,6) 

'"t?(£g!S;S^^^^^^^^ i^32(/i, 6)021 (M, 6)&32(6, + &32(^, e2)c3l(/i, ei)c23(6. ^l)} 

apply EQ.|[20ll-ll.3.2j 

I bi2 (^,$1 )C31 (^,6)^32 (^,6)^21 (gl ,6)^21 (6 ,6) 

632(6,6) 

apply Eq.|[T3t-l3.1,2> 

^32 (/^7 ^1) 632 (At,6)c2i (6 ,6)"2 (6 i6)c3i (/J, 6 )c2i (6 ,6 ) i fci2 (At,6 )c3i (^,6)632 (m,6 )fe2i (6 ,6 )c2i (6 ,6 ) 
632(^2,^1) 

apply Ea. lp^ -I3.1.2> 



fel2(/i,6; C3l(M,6)fe32(M,6) 



{ C2I (6 , ^1 ) «2 (6 , 6 ) C2I (6 , ^3 ) + &12 (6 , 6 ) C2I (6 , ^3 ) &21 (6 , ^3 ) } 



apply Y-B (13,5) 

C31 (^,6)632 (^,6)632 (m,6 )«2 (6 ,g)C2i (6 ,6)q2 (gi ,6) 
632(6,6)632(6,6) 



58 





{C32(6, 6)c2l(6, ^3)&32(6, ^3) + ^32(6, 6)c3l(6, ^3)C23(6, 6)} 
^ 



apply Eq.|[20t-ll.3.2l 

apply Eq.tT2t-f3.1.2> 



fel3(M,6)c3l(M,6)fe32(M.g3)&32(^l; ^2)^31 (Cl) Cs) £21(6,6) 
f'32(Cl,€3)''32{6.6) 




{C3l(/U, 6)a3(/W, 6)c3l(^l, 6) + &13(/U, 6)c3l(/U, 6)^31 (6, 6)} 
^ ^ 
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